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Abstract 


Year-round observations of the physical snow and ice properties and processes that govern 
the ice pack evolution and its interaction with the atmosphere and the ocean were conducted 
during the Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAIC) 
expedition of the research vessel Polarstern in the Arctic Ocean from October 2019 to 
September 2020. This work was embedded into the interdisciplinary design of the five MOSAIC 
teams, studying the atmosphere, the sea ice, the ocean, the ecosystem and biogeochemical 
processes. The overall aim of the snow and sea ice observations during MOSAiC was to 
characterize the physical properties of the snow and ice cover comprehensively in the central 
Arctic over an entire annual cycle. This objective was achieved by detailed observations of 
physical properties, and of energy and mass balance of snow and ice. By studying snow and 
sea ice dynamics over nested spatial scales from centimeters to tens of kilometers, the 
variability across scales can be considered. On-ice observations of in-situ and remote sensing 
properties of the different surface types over all seasons will help to improve numerical 
process and climate models, and to establish and validate novel satellite remote sensing 
methods; the linkages to accompanying airborne measurements, satellite observations, and 
results of numerical models are discussed. We found large spatial variabilities of snow 
metamorphism and thermal regimes impacting sea ice growth. We conclude that the highly 
variable snow cover needs to be considered in more detail (in observations, remote sensing 
and models) to better understand snow-related feedback processes. The ice pack revealed 
rapid transformations and motions along the drift in all seasons. The number of coupled ice- 
ocean interface processes observed in detail are expected to guide upcoming research with 
respect to the changing Arctic sea ice. 
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1 Gaps of knowledge 

Sea ice and its snow cover make the Arctic Ocean distinct from most other oceans; they 
control energy transfer and important interactions between the atmosphere and the ocean. 
Therefore, they require a realistic representation in Arctic climate models, which are needed 
to improve our understanding of causes and consequences of the currently observed dramatic 
changes of the Arctic climate system (Thoman et al., 2020). Although model projections are 
continuously improved and becoming more consistent in forecasting the further decline of 
the Arctic sea ice cover, they still differ considerably on the rate of the decline (AMAP, 2017). 


After the pioneering drift of Fridtjof Nansen with his vessel Fram (1893-1896) (Nansen, 1897), 
numerous expeditions were performed to study the properties and interactions of Arctic sea 
ice and its snow cover with the atmosphere and the ocean. Most of these studies were 
internationally coordinated activities using icebreakers. Furthermore, 41 Russian North Pole 
drifting ice camps were established between 1937 and 2015 (Frolov et al., 2005). Only a few 
studies were conducted observing the evolution of sea ice and snow cover during different 
seasons or over a full annual cycle. During the International Geophysical Year (1957-1958), 
sea ice heat and mass balance measurements were made at the Ice Station Alpha in the 
Beaufort Sea (Untersteiner, 1961). The Arctic Ice Dynamics Experiment (AIDJEX) focused on 
understanding ice mechanics and dynamics during the year-long drifting ice camp in the 
Beaufort Sea from 1975 to 1976 (Untersteiner et al., 2007). The international Surface Heat 
Budget of the Arctic Ocean (SHEBA) field campaign was a year-long drift of the Canadian Coast 
Guard research icebreaker Des Groseilliers in the Beaufort and Chukchi Seas from 1997 to 
1998. SHEBA focused on understanding the ice albedo and cloud radiation feedback 
mechanisms to improve climate models (Perovich et al., 1999). During the International Polar 
Year 2007/2008, snow and sea ice processes were studied along the drift of the schooner Tara 
(Gascard et al., 2008; Nicolaus et al., 2010; Haas et al., 2011). Between 2002 and 2008, 
Canadian Arctic programs involved the overwintering of the Canadian Coast Guard research 
icebreaker Amundsen to study the role of sea ice (Fortier and Cochran, 2008; Barber et al., 
2010). The multi-disciplinary Norwegian young ICE (N-ICE2015) expedition using the 
Norwegian research vessel Lance was carried out from January to June 2015. N-ICE2015 
focused on the transition processes from multi-year to younger and thinner sea ice in the 
Arctic Ocean. It spent nearly 6 months studying atmosphere-ice-ocean-ecosystem 
interactions in an ice pack dominated by relatively thin (< 1.5 m) first year ice (FYI) and second 
year ice (SYI) (Granskog et al., 2018). 


Beyond these crewed programs, autonomous drifting buoys, such as those coordinated by the 
International Arctic Buoy Programme (Rigor and Ortmeyer, 2003), observe annually the 
atmospheric, snow, sea ice, and ocean properties of the Arctic sea ice. Sea ice observations 
and data from airborne campaigns build another important source of Arctic snow and sea ice 
data (e.g., Haas et al., 2010; Kwok et al., 2017; Lawrence et al., 2018; MacGregor et al., 2021), 
supplemented by sporadic measurements collected by adventurers (Gerland and Haas, 2011). 
Numerous satellite observations provide an Arctic-wide perspective of sea ice properties 
including extent, concentration, thickness, type, drift, snow cover, and others. Passive 
microwave satellites have provided a sea ice climatology over more than 40 years, but they 
lack details for many ice-related processes (Kern et al., 2019). To establish new satellite 
products, coordinated ground truth observations are required (Gerland et al., 2019; Konig et 
al., 2019). 
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As a result of past expeditions, the understanding of the role of snow and sea ice in the Arctic 
and its influence on global climate has improved in recent decades. However, there are still 
significant knowledge gaps (Webster et al., 2018; Gerland et al., 2019) in quantifying the 
interactions and feedback mechanisms between the physical, ecological, and biogeochemical 
processes. Furthermore, an accurate representation of sea ice effects on biology and 
biogeochemistry needs to be better understood. To contribute to solving these problems, the 
work of the Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAIC) 
snow and sea ice team (in short ICE team) aims to achieve five overarching goals: 
e Characterize the properties of snow and ice cover and understand the processes that 
govern these properties, 
e Determine the snow and sea ice mass and freshwater balances, 
e Quantify the partitioning of solar radiation between the snow, the sea ice, and the 
ocean, 
e Describe the spatial variability and temporal evolution of the snow and ice cover, and 
e Integrate snow and sea ice measurements with the atmosphere, ocean, and biosphere 
of the coupled Arctic system. 


To achieve these main goals, a coordinated, integrated, and interdisciplinary approach was 
implemented as part of the MOSAIC expedition. Section 2 introduces the MOSAIC expedition 
and explains the concept of the snow and sea ice research program. The work program and a 
synopsis of the resulting data sets are presented in Section 3. In Section 4, the methodological 
and conceptual advances, and exemplary results from two case studies during one week of 
winter and one week of summer observations are shown. They result in a general description 
of the seasonal changes of the ice pack and the associated challenges of consistent and 
integrated observations of the coupled system components. Finally, linkages to the different 
Arctic sub-systems, numerical modelling, and satellite observations are discussed in Section 
5. Section 6 provides the major conclusions of this paper. 
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2 Methods and approach 

The MOSAIC field measurements started with the departure of the German research ice 
breaker Polarstern (Alfred-Wegener-Institut Helmholtz-Zentrum fir Polar- und 
Meeresforschung, 2017) from Troms@g, Norway, on September 20, 2019, and ended with the 
arrival of the vessel on October 12, 2020, in Bremerhaven, Germany. Key dates are 
summarized in Table 1 and a map is shown in Figure 1. The observational year was divided 
into 5 legs: Leg 1 included the set-up of the first central observatory (CO1, ice camp plus 
installations on Polarstern) and was supported by the Russian research ice breaker Akademik 
Fedorov. Akademik Fedorov performed most of the deployment of the distributed network 
(DN). The DN consisted of an hierarchy of autonomous systems deployed on ice floes 
surrounding the CO at 3-40 km distance to measure atmospheric and oceanic lateral gradients 
and increase sampling the highly heterogeneous ice pack properties. The network consisted 
of more than 100 simple position nodes, 8 nodes of medium instrumentation (M-sites) and 3 
nodes of large sensor suites (L-sites). A more detailed description of the DN is under 
development, led by B. Rabe. The winter Leg 2 and spring Leg 3 continued the work on CO1, 
before Polarstern had to end the (manned) Drift1 and leave the floe, for logistical reasons, on 
May 16, 2020. The autonomous stations on CO1 and in DN1 continued the (unmanned) drift 
and recorded data while the vessel performed logistical operations in Svalbard. Afterwards, 
Polarstern returned to the original ice floe, but at a different location some hundred meters 
away. Leg 4 continued the drift with the new CO2 over summer until the decay of the floe in 
Fram Strait (between Greenland and Svalbard; Figure 1) on 31 July 2020. This ended the Drift2 
and only a few autonomous devices continued the drift of CO1, CO2, and DN1. After the final 
rotation of personnel, Polarstern travelled back into the ice and started Drift3 with the set-up 
of CO3 and DN2 on August 21, 2020, near the North Pole. Drift3 ended on September 20, 
2020, when the vessel started the return voyage. A publication with details on the expedition 
in general and its logistics is under development led by M. Rex. 


Figure 1: Drift tracks of the central observatories (CO) of MOSAIC in 2019-2020. 
Table 1: Key dates of the MOSAIC expedition. 


To achieve the main goals of MOSAiC, a coordinated, integrated, and interdisciplinary 
approach was implemented, as a collaboration of five scientific teams: the atmosphere 
(ATMOS) team, the snow and sea ice (ICE, this manuscript) team, the oceanography (OCEAN) 
team, the ecosystem (ECO) team, and the biogeochemistry (BGC) team. The work on physical 
snow and sea ice processes and properties was organized in 13 tasks, which are described in 
the sub-sections of this chapter (see also Text $1 with Figure $1). All tasks comprise elements 
of observations and numerical modelling, while the balance between these elements differs 
between the tasks. Additional details on some methods are given in Text S2, including a list of 
short names and acronyms in Table S82: 

1. General snow and ice observations 
Snow measurements and sampling 
Ice coring 
Mass balance observations 
Terrestrial laser scanning 
Optical measurements 
Melt pond observations 
Sea ice dynamics 
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9. Ridge observations 

10. Transect measurements 

11. Remotely operated vehicle (ROV) observations 
12. Helicopter observations 

13. On-ice remote sensing 


The snow and sea ice measurements were distributed over the entire CO (Figure 2). A detailed 
description of all three COs with additional maps and details on installations is under 
preparation by the project team led by M. Nicolaus. While the set-up of the CO changed 
between legs due to ice dynamics and ship relocation, its general layout is exemplified in the 
map in Figure 2, which shows the CO1 ice camp shortly before the temporary departure of 
Polarstern at the end of Leg 3. 


First results with a focus on snow and sea ice conditions have been published, e.g., on the 
history of the ice drift prior to the MOSAIC expedition (Krumpen et al., 2020), on the ice 
conditions from satellite observation along the drift track in previous years (Krumpen et al., 
2021), on the representativeness of the selected ice floe (Belter et al., 2021), on the linkages 
to the general atmospheric circulation (Dethloff et al., 2021), on sea ice remote sensing 
methods (Munoz-Martin et al., 2020; Stroeve et al., 2020; Semmling et al., 2021), and on the 
platelet ice accumulation in winter (Katlein et al., 2020). 


Figure 2: Main sites and installations in the central observatory at the end of Leg 3 (CO1). 


2.1 General snow and ice observations 

The general snow and ice observations and specific activities were performed during all legs 
to document the general snow and ice conditions and to support data interpretation 
afterwards. Mapping activities of discrete measurements, installations, and operations on the 
ice were also organized in this task. Its activities included: 

(1) Acquiring timelapse photos of the surroundings of the vessel with a 360-degree panorama 
(Panomax) camera mounted above the crow’s nest. Each photo has a resolution of 
15,680x2,048 pixels and photos were taken at 20-minute intervals (72 photos per day). 

(2) Recording standard sea ice observations daily (on transits hourly) from the bridge following 
the protocol of the Arctic Shipbourne Sea Ice Standardization Tool (Hutchings and Faber, 
2018). The aim of these observations was to describe the surrounding conditions over a 10- 
minute period within a radius of 1.5 nautical miles around the vessel, including the polar night. 
(3) Documenting the on-ice work by establishing a local floe-based x, y coordinate system to 
provide the relative locations of the various sites (FloeNavi system), because measurements 
on sea ice have the particular challenge that the ice floes are in constant movement. In 
addition, geographic positions were recorded to allow reference to any other measurements, 
in particular airborne and satellites, that are not referenced to the floe. 

(4) Acquiring visible and thermal photography from helicopters and drones to provide 
information about sea ice surface conditions, to support mapping of the main floe, and to 
generate digital elevation models of the surface. Visual floe maps and digital elevation models 
were created by stitching photographs from single flights. The final images were geo- 
referenced to the local x-y coordinate system. 

(5) Upgrading, before the start of the drift expedition, the Electronic Charting System on board 
Polarstern, for receiving and displaying various satellite and weather data in near real time. 
The effort was focused on the charting of weather forecast products provided by the German 
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Weather Service. In addition, model-based ice drift forecasts were included, as well as a 
combination of various new satellite products. 

(6) Documenting the on-ice work almost daily on board and synchronizing with land as the 
‘MOSAIC LogBook’ to inform the project consortium on details of the ongoing work. 


2.2 Snow measurements and sampling 

The overarching goal of the snow measurements was to characterize the spatial and temporal 
variability of snow over different sea ice types (e.g., FYI, SYl and ridges) during an entire annual 
cycle. During the melt season, techniques used to investigate the snow were applied to the 
uppermost layer of the melting ice cover, the so-called “surface scattering layer” (SSL). The 
snow measurements during MOSAIC were conducted with the aim to improve snow models, 
and subsequently to improve climate models and remote sensing retrieval algorithms, and in 
support of all MOSAIC teams. Snow is one of the most insulative naturally occurring materials 
on Earth (Webster et al., 2018), and one of the most efficient reflectors of incident solar 
radiation (Warren et al., 1983). Snow on sea ice therefore acts as an effective barrier between 
the atmosphere and the ocean, preventing both heating of the upper ocean in summer and 
cooling of the ocean in winter. While the bulk physical properties of snow on sea ice (e.g., 
snow depth, snow density, snow thermal conductivity, snow reflectivity, albedo) and the 
resulting stratigraphy have been studied previously, data on the spatial heterogeneity of the 
optical and thermal properties of snow on sea ice and the relationship between these 
properties and their governing processes are limited. Snow measurements during MOSAIC 
ranged from microscale (micrometer) to local scale (approximately 1000 m) to address these 
properties across scales. The following novel instruments were used during MOSAIC to 
measure physical properties of snow: 

- The Snow Micro Pen (SMP), a penetrometer (Proksch et al., 2015) for fast retrieval of 
the vertical profile of snow stratigraphy, snow density, and specific surface area, with 
a high vertical resolution (better than 1 mm). The SMP measures vertical profiles of 
snow layers at one location. By repeating samples over a large area, micro-, macro- 
and regional scale properties can be linked. SMP was used on snow profiles and 
transects. Overall, 6,959 SMP profiles were measured. 

- AnX-ray microtomograph (micro-CT, Scanco Medical micro-CT 9; Calonne et al., 2014), 
which was installed on board Polarstern and allowed measurements of snow 
microstructure on quasi in-situ samples. We used the micro-CT to measure the full 
snow profiles throughout the year, as well as on the SSL, lead and pond ice, and sea 
ice cores. Overall, 614 samples were analyzed with the micro-CT. 

- Near Infrared (NIR) photography (Matzl and Schneebeli, 2006) to characterize the 
specific surface area of snow profiles as well as the snow and SSL surfaces. During 
MOSAIC, photographs at 850 nm and 940 nm were taken for each profile/surface 
scene. Overall, some 600 profile/surface scenes were photographed using NIR. 

- Structure from motion photograph sets to derive elevation data. Based on these data, 
high spatial resolution digital elevation models were created at each snow pit location, 
and sometimes along transects, of the surface and of the snow-ice interface. These 
“mini digital elevation models” will be used to determine small scale (millimeter) 
surface roughness that is relevant for optical modelling and satellite retrievals, and link 
data to coarser-scale digital elevation models obtained by terrestrial and airborne laser 
scanning. Overall, some 300 structure from motion sets were taken during MOSAIC. 

Snow observations were collected mainly in specific areas assigned to snow measurements. 
Those areas included the prevalent types of sea ice, namely FYI (e.g., Snow1), SYI (e.g., Snow2), 
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ridges and new ice in leads, as highlighted by yellow patches in Figure 2. Additional 
measurements were performed across the COs to connect physical snow properties to 
relevant processes that were studied by other teams. The team members conducted about 
600 snow pit measurements of different types: 
“Type C” snow pits characterized the vertical profile of the snowpack and consisted of 
SMP measurements only. 

- “Type B” snow pits included density, snow water equivalent (SWE), temperature, and 
salinity measurements in addition to observations of Type C. 

- “Type A” snow pits included the full suite of snow measurements, adding snow 
sampling (see below), micro-CT scans, SMP transect (to either side of pit), NIR 
photography, structure from motion pictures, and snow sampling for a wide range of 
constituents, in addition to observations of Type B and C. 

For practical reasons, this perfect scheme of snow pit types was not always followed in every 
pit. We performed a comprehensive snow sampling program, which was centrally organized 
and coordinated across all scientific teams. Snow samples co-located with the physical 
measurements were taken for analysis of salinity, water and sulphur stable isotopes, major 
ions, halogens and halocarbons, beryllium, black carbon, marine sugars, organic acids, 
microplastics and ice-nucleating particles. Snow sampling was combined with snow pit 
measurements whereby the frequency and number depended on the purpose of the snow 
sample. This approach ensured that all samples from each sampling event were taken at the 
same site(s) and in a consistent and reproducible manner facilitating the co-location with the 
physical properties of the snow. In addition, dielectric permittivity measurements were made 
at 50 MHz to invert snow wetness during melting conditions. 


Figure 3 shows the temporal distribution of snow pits. Distributed observations of physical 
snow properties on transects, and of snow depth (Section 2.10), and observations of snow 
surface position with terrestrial and airborne laser scanner (TLS Section 2.5 and ALS Section 
2.12) will be used to upscale these point observations. Collectively, these measurements 
provide a comprehensive characterization of the snow cover, enabling synthesis studies on 
snow evolution, spatial variability, energy budget and heat fluxes, remote sensing algorithm 
development and towards improved parameterizations in climate models. 


Figure 3: Snow and sea ice observations during the field phase of MOSAIC. 


2.3 Ice coring 

The primary aim of the ice coring work was to capture the seasonal evolution of sea ice 
physical properties over a full year, for both level FYI and level SYI at 5-cm vertical resolution, 
and to document the differences between FYI and older ice. This work was performed at two 
main coring sites (Figure S2 with Tables S2 and S3), one on FYI and the other on SYI, which was 
the oldest ice at the start of MOSAIC. We used the approach recommended by Oggier et al. 
(2020) to section the ice cores, relative to both the ice surface and ice bottom (Section 4.2), 
to better capture the processes at the surface and bottom interfaces. The physical properties 
affect energy and matter fluxes and ecosystem processes in sea ice, and are also important to 
the remote sensing of sea ice (Section 2.13). The ice coring work captures the seasonal 
evolution of level ice thickness, temperature, salinity and density. At each coring event, 
temperature profiles were measured using a thermistor probe into holes drilled at 5 or 10 cm 
intervals within 2-3 minutes of core extraction. Two cores collected within 50 cm of the 
temperature core were cut in sections and stored in sealed containers for transport to the 
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ship, after noting total length and fractures. Onboard, bulk salinity was measured on melted 
ice samples, and density was measured using the hydrostatic density method (Pustogvar and 
Kulyakhtin, 2016). The brine volume fraction can be derived from these measurements. Two 
additional cores were collected for textural and microstructure analysis: one for thin sections 
and asecond for micro-CT scans (Section 2.2). Bottom and surface sections (in most cases the 
lower and uppermost 10 cm of ice cores) were imaged on board using the micro-CT on 
Polarstern. The remaining ice cores were archived for more detailed studies of texture and 
microstructure at home laboratories after the expedition. In the sunlit period, ice cores were 
also collected to measure the optical properties at the home lab (Section 2.6). Coring was 
coordinated with the snow task, and often the coring was associated with a simple snow pit 
(which typically included temperature, density, salinity and sampling for oxygen isotopes; 
Section 2.2) supporting the interpretation of the ice core data. Coring was conducted in close 
collaboration with the ECO and BGC teams. 

In addition to the time series data from the main coring sites (Figure S2), the BGC team also 
observed the temporal evolution of both FYI and SYI at a number of additional sites (primarily 
salinity and temperature) coincidently with their own work. In addition, sediment laden ice 
cores were collected for analysis jointly with the BGC team and an occasional ice core was 
collected for analysis of microplastics. These ice cores complement the observations made by 
the ICE team and provide more information on spatial variability. Coring was also done 
repeatedly in ridges (Section 2.9) and for ice mechanical properties (Section 2.8). 
Opportunistic sampling was done especially for thin young ice in leads and linked to remote 
sensing work during some events and occasionally in the DN. The collected data set captures 
the seasonal evolution of FYI, from early growth and new saline ice to mature FYI and the rapid 
desalination in summer (Section 4). While the FYI that had survived the previous summer melt 
became SYI at the beginning of the campaign, the refreezing and new ice growth at the bottom 
of SYI could be observed until its summer decay. An overview of the ice coring activities by the 
ICE team is shown in Figure 3, and a list of existing ice cores is referenced in the data 
availability section. 


2.4 Mass balance observations 

The goal of the mass balance observations was to document the amount and spatial variability 
of ice growth during winter and surface and bottom ice melt during summer. Mass balance 
observations provide insight on the relative contributions of the atmosphere and ocean to 
changes in the sea ice thickness and snow depth. Multiple sites sampling new ice, FYI, SYI, 
deformed ice, undeformed ice, and ponded ice were included to comprehensively address the 
spatial variability of sea ice mass balance. The in-situ mass balance measurements were made 
both autonomously and manually with 37 digital temperature chains, 4 sea ice mass balance 
buoys (IMB) of different types (e.g., Polashenski et al., 2011; Jackson et al., 2013), 20 Snow 
Buoys (Nicolaus et al., 2021), and 120 manual ice thickness gauge/ablation stake pairs. The 
temperature chains and IMBs recorded vertical profiles of temperature of the air through the 
snow and ice into the upper ocean, as well as (for some IMBs) the positions of surface and 
bottom interfaces using acoustic rangefinders. Analysis of temperature chain data manually 
and by algorithms are also used to determine the position of the air-snow or -ice and the ice- 
ocean interfaces (Hoppmann et al., 2015; Cheng et al., 2020). Manual measurements at the 
ablation stake sites also provided a time series of the sea ice surface and bottom position 
changes. Thicknesses measured from ice cores collection and laser scanning measurements 
contributed to the snow and sea ice mass balance measurements (Sections 2.3, 2.5 and 2.12). 
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Together, these multiple methods of observing ice surface and bottom position are used to 
calculate time series of snow depth and ice thickness, ice growth, surface melt, and bottom 
melt. In combination with transects (Section 2.10) this approach provides a comprehensive 
view of the sea ice mass balance (in space and time) and the best local (thermodynamic) and 
best spatial (dynamic) information. To complement the mass balance observations, systems 
for atmosphere and ocean measurements were co-located on the same sites. 


2.5 Terrestrial laser scanning 

The goal of the terrestrial laser scanning (TLS) was to map the snow and ice surface 
topography and how it changed in time and space throughout the seasonal cycle. Using the 
eye-safe Riegl VZ1000 TLS during winter and the non-eye-safe Riegl VZ6000 TLS during 
summer allowed large aerial scans. Analysis of these data is being used to link these changes 
with other in-situ measurements of snow and ice properties and spatially extend point 
measurements of snow and ice properties. TLS measurements were used to map an area of 
approximately 0.5 km? every two weeks with a centimeter-scale vertical resolution and 
decimeter-scale or better horizontal resolution. The same region was mapped with 42 scans 
during the drift with CO1 from October 18, 2019, to May 9, 2020, (Figure 3) and included most 
of the snow pit measurement sites, most of the mass balance observation sites as well as the 
southern loop of the transect, all of the albedo transects in spring, portions of the strain gauge 
installations, the area studied from beneath the ice at ROV sites (ROV 2.0 and 3.0 in Figure 2), 
the on-ice remote sensing instruments, and the on-ice meteorological installations at Met City 
(Sections 2.2, 2.4, 2.6, 2.11, 2.13). This co-location of the various methods will allow linked 
data processing, while the impact of surface destruction on the TLS results through other 
activities is minor given the very small fraction (< 1%) of impacted surface area. Summertime 
TLS observations were more limited due to safety and logistics but included areas of snow 
pits, albedo, ponds, and other surface property measurements. The TLS observations are used 
to generate a time series of ice dynamics, changes in atmospheric form drag, and snow 
accumulation, redistribution and melt. TLS measurements can be used to infer material 
properties throughout the mapped region by upscaling the vertical profile measurements of 
snow and ice physical properties made by other instruments. 


2.6 Optical measurements 

The objective of the “Optics” task was to quantify the partitioning of solar radiation within 
and beneath the ice cover during the sunlit season. The quantities measured included albedo, 
transmittance, and light extinction in snow, sea ice and the underlying water. A variety of 
sensors (measuring either broadband, spectral or photosynthetically active radiation) were 
used for measuring repeat-visit survey lines and grids in addition to fixed-point study sites. 
The albedo measurements included surface survey lines, multiple fixed site installations, 
unmanned aerial systems (e.g., the drones HELIX and Mavic and the Spectra Quadrocopter), 
helicopter-borne instrument platforms and the tethered balloon. Transmittance 
measurements included repeat under-ice surveys using an ROV (Section 2.11) and site-specific 
measurements using an L-arm photosynthetically active radiation sensor. In-ice extinction and 
inherent optical property measurements were carried out at specific locations using profiling 
sensor arrays frozen into the ice (light chains and light harp) and a spatially resolved 
reflectance probe lowered into auger holes. Light extinction within the water column beneath 
the ice cover was estimated using sensors mounted on the ocean CTD, and the ROV was also 
equipped with a number of optical sensors. 
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Sensors were also deployed on radiation stations to autonomously collect albedo, 
transmittance and water column radiation at fixed locations throughout the season. In 
addition to these time series, specific experiments were carried out to assess the response of 
spectral albedo to specific events or processes: 1) snow during early summer melt, 2) diurnal 
variability during both melt and freeze-up, 3) the removal and formation of a surface 
scattering layer, and 4) rain on new snow. Additionally, ice was sampled for investigation of 
the linkages of optical and structural properties to be carried out in freezer laboratories at 
home institutions (Section 2.3). 


Ongoing investigations of the spatial heterogeneity of the surface albedo at various horizontal 
scales will be aided by the drone-based spectral and broadband albedo measured along grids 
and horizontal transects at various heights above the ice surface. The vertical profiles of 
albedo recorded from drones will be crucial for determining the contributions of various 
surface features (melt ponds, ridges, flat-ice, leads) to the areal-averaged albedo. These 
measurements enable scaling from ground-based observations (meter-to-tens-of-meters 
footprint), to satellite observations (tens-to-hundreds-of-meters footprint) and large-grid 
model outputs. Moreover, high spatial resolution hyperspectral albedo measurements 
collected from drone-mounted cameras will help relate the scattering properties of each 
surface feature to the snow and ice microstructure properties observed at the surface. 


2.7 Melt pond observations 

Melt ponds play a key role in energy exchange between atmosphere and ocean in the Arctic 
summer because they substantially reduce surface albedo of sea ice. Compared to bare ice, 
melt ponds reflect less solar energy to the atmosphere and transmit more solar energy into 
the ice underneath and the ocean, thus fostering further sea ice melt (Light et al., 2008; 
Nicolaus et al., 2012). At present, a synoptic view of pond evolution, distribution and depth 
using airborne and satellite data is still limited. The main objective of the Ponds task was to 
provide field and airborne data (Section 2.12) of the physical pond characteristics (e.g., size 
and shape distributions) to improve the understanding of the spatio-temporal development 
during all stages of pond evolution. These characteristics include pond area fraction, size, 
shape, and number. To improve pre-processing of airborne and satellite data (i.e., 
atmospheric correction, correction of adjacency effects; K6nig et al., 2019), the spectral 
signature of the different surface types on the floe were measured using hand-held water 
color spectroradiometers, a hyperspectral camera, and a goniometer. Up- and downwelling 
broadband irradiance data as well as multi-spectral imagery are available from drone surveys, 
helicopter flights and satellite observations (Sections 2.6 and 2.12). To assist a proper 
reflectance calibration of water surfaces we used a remotely controlled boat equipped with 
water color spectroradiometers. Above and below water spectral measurements were used 
to assess the optical properties of pond water (attenuation coefficient, inherent optical 
properties; K6nig and Oppelt, 2020). In addition to the radiometer measurements, the field 
measurements also included pond depth as well as chlorophyll content, salinity and 
temperature of pond water. To assess the absorption coefficient of the dissolved and 
particulate matter, water samples were taken and analyzed in the Polarstern lab using a point- 
source integrating-cavity absorption meter. Micro-CT samples were processed from frozen 
melt ponds to derive optical properties. 
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2.8 Sea ice dynamics 

One of the key questions of MOSAIC is how sea ice moves and deforms. The dynamics and 
mechanics of the ice around the Polarstern and its interaction with the ship, ocean and 
atmosphere were monitored throughout the drift. We aimed to examine the spatial and 
temporal distribution of ice motion and ice stress, and to determine how they vary locally, 
regionally, and seasonally in relation to the physical properties of the ice. The CO observations 
were coordinated with satellite and airborne remote sensing (Section 2.12) monitoring ice 
drift, ice morphology and floe size. Drift and strain-rate (divergence, shear and vorticity) of 
the ice pack were measured over multiple nested scales, from hundreds of kilometers in a 
wide array of buoys coordinated with the International Arctic Buoy Program, to ocean 
synoptic and mesoscale scale in the DN measured with more than 100 ice-deployed GPS 
position sensors, to ship-mounted ice radar (Oikkonen et al., 2017) providing dense target 
tracking over 7 km around the ship every few minutes, to high frequency and centimeter-level 
position accuracy from GNSS helicopters inertial navigation system stations and robotic 
theodolite strain observations. Ice motion was also monitored at two locations with seismic 
stations, providing information on vertical motion, and ice-wave interaction (Marsan et al., 
2011). Three seismic stations, located at three points of the ice floe, recorded ice oscillations 
in three directions (horizontal and vertical) and transmitted the records to the central station 
installed on the ship. The signal recording frequency was 100 Hz. The graphical format of the 
waveform records allows to see ice events such as new fractures, compression or ridging. 
Internal ice stress was recorded with vibrating wire strain gauges (Cox and Johnson, 1983) at 
17 sites in the CO. Twelve of these gauges had 1-minute sampling, the others 10-minute 
sampling, allowing for specific breaking events to be aligned with the stress time series. To 
estimate the ice stress on the Polarstern, stress panels were placed in the ice about 20 m from 
the ship parallel to ‘Void 92’, a specific compartment in the hull of Polarstern, which was 
instrumented with strain gauges. Throughout the drift the mechanical properties of the ice 
(local strength, horizontal and vertical uniaxial strength tests) were measured at a level SY| 
site. For these measurements, an LGK hydro-complex with a probe-indenter and a press was 
used. The probe-indenter was used to measure the local strength of the ice cover at different 
levels, and the press was used to measure the ice strength for uniaxial compression (horizontal 
and vertical). Salinity, temperature, density and texture analysis of co-located ice cores allow 
characterization of the mechanical properties of ice (Section 2.3). A log of specific dynamic 
events (cracks, leads opening, ridging) was kept, Panomax camera and Ice Radar were used to 
identify the start time of each event to within 20-minute precision. The full data set allows for 
specific case studies that can correlate ice fractures and the response of ice to ocean and wind 
forcing on scales beyond tens of meters. 


2.9 Ridge observations 

Ridges are one of the least studied features in the Arctic ice pack, because of logistical 
challenges in making measurements in these sometimes-massive chaotic formations of ice 
that pile up when ice floes collide. The main aim of the ridge work was to study the 
development of the geometry and of the thermo-mechanical properties of ridges during the 
MOSAIC drift. The three main phases of ridge evolution, initial formation, winter consolidation 
and summer decay, were documented for the first time for selected ridges throughout a full 
seasonal cycle and compared to the evolution of level FYI and SYI. The main geometrical 
parameters characterizing the ridges were the thickness of the consolidated layer, the keel 
depth and the macro-porosity. These three parameters were determined by repeated manual 
drilling (Figure 3) and complemented by continuous temperature profiles obtained by IMBs 
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and thermistorstrings that were deployed vertically through specific ridges. Additionally, sea 
ice cores were taken from ridges to characterize the temporal and spatial evolution of sea ice 
temperature, salinity and density. In addition to this baseline work, additional specialized 
measurements were performed to determine bottom topography, ocean properties, currents 
and keel-associated turbulence, as well as surface characteristics by snow and ice thickness 
transects, snow surveys and laser scans of surface topography (Sections 2.2, 2.3, 2.5, 2.10 and 
2.11). The characterization of the physical properties of ridges was of particular importance 
to coordinated ecological studies of sea ice ridges performed by the ECO team. 


2.10 Transect measurements 

Making systematic and representative observations of physical sea ice parameters, their 
spatial variability and temporal evolution in the CO throughout the full MOSAIC drift was the 
objective of the transect task. The sampling strategy consisted of repeat surveys of sea ice and 
snow thicknesses throughout all MOSAIC legs, similar to observations carried out from May 
to December during the Tara expedition in 2007 (Haas et al., 2011). These baseline 
observations were complemented with seasonal activities, specifically measurements of 
surface albedo and melt pond properties during summer and autumn. In addition, cross-task 
activities with joint measurements by mobile on-ice remote-sensing sensors (Section 2.13; L- 
band microwave radiometer and Ku/Ka-band radar) along the transect lines or transect 
observations at specific ridges (Section 2.9) were carried out to provide input data for 
algorithm development and temporal and spatial evolution at ridges. 

The transect measurements were carried out by repeating dedicated loops typically 1 km or 
longer, except for Leg 5 which was shorter. The baseline observations included a broadband 
electromagnetic induction sensor (Geophex GEM-2, Hunkeler et al., 2016) for the 
measurement of the total thickness of the sea ice and snow/SSL layers and an automated 
snow/SSL depth probe (Snow Hydro MagnaProbe; Sturm and Holmgren, 1999) which in 
combination yield point measurements of snow/SSL and sea ice thicknesses. To accommodate 
measurements of snow/SSL and melt pond depths when both were present, the snow disc of 
the MagnaProbe was outfitted with floatation on its upper side to also measure pond depth. 
During winter, only limited FYl and mixed type transects were performed, as no observations 
on FYI were possible in early winter. These transects were not revisited on later legs as they 
became inaccessible due to deformation. During summer and autumn, additional transects 
were conducted parallel to albedo lines for relating the properties of snow, SSL, sea ice, and 
melt ponds to the spectral and broadband albedo measurements (Sections 2.6 and 2.7). 

The two goals for the transect surveys were to achieve representative coverage of sea ice 
conditions in the CO and continuity of observations throughout the drift. However, due to the 
relocation of Polarstern after Legs 3 and 4 (Table 1), new transect loops were established in 
addition to the original transect loops. In the beginning of the MOSAIC expedition, the two 
original loops were laid out, termed “northern” and “southern” loop. The rationale behind the 
two loops was to survey the two dominant ice types in the CO: the northern loop on the older, 
thicker, and more deformed part of the “MOSAIC floe”; and the southern loop on the younger 
and thinner part that had widespread refrozen melt ponds. While these two loops were 
traversed continuously from October 2019 to early May 2020, the southern loop was often 
inaccessible after sea ice deformation events. The summer transect (CO2) was established in 
June 2020 on similar ice conditions as the prior transect loops and with small (approximately 
300 m) overlap with the northern transect loop. The transect loop of Leg 5 (CO3) was fully 
independent of earlier transects due to the change in regional focus. A full overview of the 
dates, location, and acquisitions of all transect observations is referenced in the data 
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availability section. In addition, the dates of all total ice plus snow thickness observations 
including data acquisitions outside the transect task are illustrated in Figure 3. 


2.11 Remotely operated vehicle (ROV) observations 

The main goal of the work with the ROV was to obtain a better understanding of the seasonal 
cycle with respect to (1) ice draft and bottom topography, (2) the light field beneath the ice, 
(3) the bio-physical properties of the ice and uppermost ocean, and (4) the organisms living in 
and under the sea ice. All measurements focused on the spatial variability under sea ice and 
across different surface features and ice types. To accomplish these objectives, the ROV was 
operated one to three times per week amounting to a total of 83 dives. The ROV was equipped 
with a comprehensive suite of bio-physical sensors, several cameras, and a mechanism to haul 
sampling nets under the ice (Katlein et al., 2017; Wollenburg et al., 2020). During summer, 
additional dives were performed with an underwater hyperspectral imager (UHI) to survey 
the 2D and 3D radiance under sea ice. In addition to the direct data sets from the ROV, the 
work was complemented by data sets from other ICE and ECO team tasks to allow the up- 
scaling of point observations, observing links to ecosystem parameters, and estimating energy 
budgets combined with other observations. 

The ROV was operated within an approximate 300 m radius around the ROV sites (Figures 2 
and S3). Over the year, the ROV site was moved several times due to the dynamic icescape. 
Specific dive missions were defined and executed as part of the weekly schedule. Over most 
of the year, weekly multi-beam sonar surveys and net hauls were carried out. Multi-beam 
sonar surveys mapped the under-ice topography in a grid covering an area as large as possible 
at a depth of about 20 m (Section 4.1). Net hauls were conducted along linear or triangular 
transects at different depths for dive times of approximately 15 minutes. Once sufficient 
sunlight returned (mid-March), comprehensive optical dives were carried out on the ROV days 
(including downwelling radiation), in particular for grids and transects under FYI, SYI, leads 
and ridges. This effort also included stable measurements at marker positions, vertical profiles 
and measurements of upwelling radiation. Additionally, specific dive missions were added to 
locate lost devices, to document under-ice installations and to deploy and recover sediment 
traps. 


2.12 Helicopter observations 

The main objectives of the helicopter-borne sea ice surveys were, first, to document the 
temporal evolution of snow and sea ice cover on the CO floe and within the DN, from sea ice 
formation in autumn to sea ice melt in summer and, second, to document the spatial 
distribution and spatial variability of snow and sea ice properties within the DN and beyond. 
Key parameters were sea ice thickness distribution, surface topography, freeboard and floe 
size distribution, surface temperature and surface albedo, areal fraction of different sea ice 
types, and melt ponds as well as additional melt pond characteristics like pond-size 
distribution and bathymetry (Section 2.7). 

To reach these goals, three different sensor packages were operated independently on 
different flights: (1) an airborne laser scanner (ALS, Riegl VQ-580) , two RGB cameras (Canon 
EOS 1D Mark III) with wide-angle and fisheye lens, in winter an IR camera (VarioCAM HD head 
680, InfraTec) or in summer a VIS/NIR hyperspectral camera (Specim AisaEAGLE, 400 nm to 
970 nm), a radiation thermometer (Heitronics KT19) and an upward and downward looking 
pyranometer (Kipp&Zonen CMP22), (2) EM-Bird for measuring the combined thickness 
distribution of sea ice and snow using airborne electromagnetic induction sounding (Haas et 
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al., 2009), and (3) the towed system HELiIPOD, which was operated by the BGC team for 
measurements of trace gases, energy and momentum fluxes, aerosol and sea ice surface 
properties. The laser scanner and camera configuration was operated during all legs on a total 
of 68 flights, the EM-Bird was operated during all legs, except in the midwinter darkness, on 
23 flights, and the HELIPOD was operated in spring and summer on five flights. 

There were five main mission types to investigate the temporal evolution and spatial 
variability of snow and sea ice properties: (1) local floe grid survey covering the CO floe and 
surrounding floes, (2) regional survey consisting of lines between the three L-site buoys and 
between these buoys and Polarstern, (3) butterfly survey consisting of triangles between 
Polarstern and four selected buoys of the DN to observe the relation between ice deformation 
and thickness change, and (4) mixed pattern of horizontal lines and vertical profiles mostly 
across the CO. Missions (1) and (4) were flown with the ALS and camera sensor package as 
well as with the EM-Bird and HELiPOD, mission (2) with the ALS and camera sensor package 
only, and mission (3) with the EM-Bird only. The maximum airtime of individual flights was 
about 2 hours. 


2.13 On-ice remote sensing 

Almost all Arctic-wide sea ice climate time series are based on microwave satellite 
observations (IPCC, 2013), which now span more than 40 years. However, the retrieved 
geophysical parameters like ice concentration, ice type, ice thickness, melt onset, or snow 
depth are not measured directly by the microwave radiometers and radars but inferred from 
the interaction of microwave radiation with sea ice and snow as well as seawater. Depending 
on frequency and physical properties (temperature, salinity, snow grain size and geometry, 
porosity, etc.) of the snow and sea ice (including melt ponds), microwave radiation can 
penetrate into the ice, get scattered, absorbed and emitted. Information of the microwave 
radiative transfer in snow and sea ice is needed to derive the observational quantities from 
satellite measurements. 


The main objective of the ice remote sensing (RS) task is to develop new methods to retrieve 
improved sea ice and snow parameters from satellites. This goal relates both to reducing and 
quantifying the uncertainties of existing satellite retrievals (e.g., ice concentration, sea ice 
thickness and ice type) and to developing new methods for current and future satellite 
missions (e.g., snow depth for CRISTAL and CIMR, surface properties and ice types from 
Radarsat Constellation Mission and TerraSAR-X, melt pond depth using hyperspectral airborne 
data and Sentinel-2). Most measurements described here were in the microwave domain 
from 0.5 to 89 GHz using passive microwave radiometers and active radars. Here, the most 
pressing issue is to better understand how the microwave emissivity and backscatter depend 
on physical snow and ice properties and environmental conditions. Thus, remote sensing 
instruments were co-deployed with snow and ice measurements (Figures 2 and S4). In turn, 
radiative transfer modelling of snow and ice will be improved with the availability of the RS 
measurements. Further observations were carried out in the visual and IR domain to support 
the microwave measurements and to measure ice, snow and melt pond optical properties as 
well as surface temperature distribution and evolution. MOSAIC offered the unique 
opportunity to observe a full seasonal cycle and co-locate measurements at different 
frequencies and polarizations. RS measurements were taken at a central remote sensing site 
on the ice floe during all five legs, from the ship, and along transects (Section 2.10). An 
overview of the remote sensing site is given in Figure S4. A detailed overview of satellite 
related work is under development led by G. Spreen. Sea ice conditions along the drift 
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trajectory based on satellite observations and comparison to previous years are described in 
Krumpen et al. (2021). 


Table 2 provides a short overview of the 14 RS instruments that were operated during 
MOSAIC, with their basic instrument parameters. In addition, examples of corresponding 
satellites or sensors are given. The three main instrument categories were microwave 
radiometers, radars, and other sensors like GNSS-R, IR, and cameras. These remote sensing 
measurements on the ice correspond to a wide range of satellite sensors, which are also listed 
in the table. Due to ice dynamics and instrument problems not all instruments were 
operational all the time (see Figure 3b for observation periods). The primary measured 
parameters are emitted infrared and microwave radiation (brightness temperature from 
radiometers), backscattered (radars) and reflected (GNSS-R, nadir Ku/Ka radar) microwave 
radiation, and spectral reflectances (hyperspectral camera). 


Include Table 2 here 
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3 Synergy in the observational program and available data sets 

The work on snow and Sea ice processes and properties required an integrated approach to 
carry out the observational program along the tasks outlined in Section 2. Performing 
consistent measurements throughout a full year was a major challenge and one of the main 
differences between MOSAIC and many previous snow and sea ice field studies. To guarantee 
consistent high-quality data collection between all legs, a particular preparation phase was 
performed (Text $1). This section describes the key processes, which were studied by the ICE 
team across the different tasks (Section 2). These processes are illustrated in Figure 4. While 
most processes were observed continuously over the full annual cycle, seasonal processes 
were only studied during parts of the drift (e.g., solar short-wave radiation interactions). The 
work covered multiple ice types and the highly variable snow distribution: new ice, FYI, and 
SYI. We worked on level, deformed and ridged sea ice as well as in the marginal ice zone, leads 
and melt ponds. In addition, we describe the measured parameters and the resulting data 
sets, which will be used to improve and develop parameterizations for a better understanding 
of the coupled Arctic sea ice system. The key parameters are summarized in Table 3 and their 
temporal coverage is shown in Figure 3. Methods and tasks are restricted to the work in the 
COs. 


Figure 4: Schematics of the snow and sea ice processes studied during MOSAIC. 
Include Table 3 here 


Sea ice and snow mass balance 

New ice formation in open water (mainly in leads) was observed through the entire autumn, 
winter and spring, including observations of surface, bottom, lateral and internal freezing. In 
October 2019, most measurements were performed on SYI that had survived the previous 
summer. However, parts of the ice camp, e.g., the coring activities, were also established on 
younger sea ice that involved FYI. Thus, the evolution of both FYI and SYI during the initial 
freezing was observed and followed throughout the campaign. Furthermore, MOSAIC was 
able to observe the initial accumulation of snow in autumn. Snow accumulation, distribution 
and re-distribution were observed through the entire year with a particular focus on 
comparisons of snow processes between level FYI, level SYI, and deformed (ridged) ice in the 
CO, and how these processes affect the growth of different ice types. 


The novel observation of platelet ice under winter sea ice in the Arctic showed that platelet 
ice contributed to bottom ice growth during winter (Katlein et al., 2020). Over time, the ice 
pack became strongly deformed. After melt onset, the melt and decay of the snow and sea 
ice were observed, and contribution from surface, bottom, lateral and internal melt processes 
were quantified. The decay of the original MOSAIC floe (CO1 and CO2) on July 31, 2020, 
showed strong linkages between dynamic and thermodynamic processes controlling the 
decay of the ice pack in the marginal ice zone. 


In addition to the general time series of each ice type, the comparison of sea ice and snow 
mass balance for level and deformed ice was an important element along the transects. The 
thickness, consolidation and decay of pressure ridges were examined by a combination of 
drilling, coring and analysis of IMB and thermistor string data. The macro-porosity of ridges 
was estimated from drillings while the keel depth was investigated by drilling, ROV multi-beam 
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mapping and IMB/thermistor data. Almost all tasks contributed in some way to the sea ice 
and snow mass balance work (Table 3). 


Ice and snow physical properties 

The physical properties of sea ice were studied based on ice cores and samples on the micro- 
scale, but also through surface and under-ice transects as well as larger scale airborne and 
dynamical studies (Table 3). The ice coring work of the ICE team focused on the time series of 
the main coring sites to capture the evolution of FYI, SYl and deformed ice. Additional 
information on the spatial variability is available through the work related to the mechanical 
properties, the remote sensing work and in particular through the BGC team, which surveyed 
additional sites across the COs and DN1. Optical properties were mostly studied in relation to 
energy budget studies (see below). 


Snow on sea ice is the thermal and structural interface between sea ice and the atmosphere 
(Sturm and Massom, 2017). Thus, even though major parts of the snow work were performed 
by the ICE team, many linkages exist to the ATMOS team, and comprehensive parts of the 
snow work were carried out in strong collaboration with the ATMOS team. Snow precipitation, 
rain fall, and snow re-distribution were studied with common set-ups by both teams. Snow 
formation, snow metamorphism and changes in the snowpack were studied in detail from 
micrometer scales to floe scales. Snow cover is the most (temporally and spatially) variable 
layer dominating the thermal and optical properties of the sea ice system. Due to the highly 
metamorphic nature of snow, its microstructure changed at a fast pace. Snow processes 
during MOSAIC were highly sensitive to quickly changing air temperatures, which led to strong 
temperature gradients in the snow, particularly during Legs 3 (spring) and 5 (autumn), and 
resulted in a high degree of snow metamorphism. The strong temperature gradients affected 
the thermodynamics, growth rate and microstructure of the underlying sea ice, and 
consequently strongly affected the temperature, salinity and permeability of sea ice. 
Ultimately, the physical snow properties and processes build the link to light and nutrient 
availability in the sea ice and in the upper ocean, and significantly affect the polar marine 
ecosystem. 


Energy budget, optical properties, and melt ponds 

Incident solar shortwave radiation is partitioned into reflected, absorbed, and transmitted 
components during the sunlit season. This split determines the surface shortwave radiative 
budget, with effects on mass balance and physical properties of snow and sea ice, and 
contributes to ocean heat (Figure 4). The optical properties of the ice cover were investigated 
by a variety of methods. This variety was needed to constrain radiation budgets, provide some 
spatial and temporal overlap, achieve redundancy to help alleviate data gaps where 
instruments failed, and to provide insight into physics that are driven by or drive processes 
not well constrained by a single set of measurements. An example would be the utility of 
surveying both albedo and transmittance simultaneously to understand the full partitioning. 
While the work of the ICE team focused on studies of this partitioning within the snow and 
ice, the ATMOS team observed radiative transfer in the atmosphere as well as the full surface 
energy budget including radiative, turbulent, and conductive heat fluxes. The OCEAN team 
observed radiative heating and turbulent mixing in the ocean, using a combination of near-ice 
direct heat, salt and momentum ocean flux measurements and profiling CTD and 
microstructure dissipation profiles. 
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Conductive heat flux between atmosphere and ocean through the snow and sea ice pack was 
studied mostly based on thermistor chains (Section 2.4), but also as residuals from surface 
and basal mass balance calculations. These heat fluxes are crucial for quantifying melt and 
freeze processes at both interfaces and thus are key elements to close the atmospheric and 
oceanographic energy budgets. Heat fluxes through snow and sea ice strongly determine sea 
ice mass balance. 


Fixed stations measured the temporal variability at single locations. Surveys measured both 
spatial and temporal variability. Fixed installations included automated radiation stations, the 
UHI, light chains hanging through the ice into the ocean below the ice, and a light harp, which 
is a novel in-ice light profile sensor that was co-located with a salt harp system. Optical 
property survey work included measurements of surface and drone-based albedo and 
through-ice light transmittance recorded on the ROV. Other measurements that did not use a 
repeat-survey sampling strategy, but also were not tied to a fixed field of view on the ice, 
included helicopter and HELiPOD albedo, in-ice inherent optical property probe, and ice cores 
taken for dedicated laboratory assessment of inherent optical property profiles (Table 3). 


The surfaces of snow and sea ice, which comprise the dominant properties for remote sensing 
and large-scale linkages, changed seasonally. While our work during winter was strongly 
related to changes in snow accumulation, ridging, and surface roughness, surface melt 
processes were examined in detail over summer. Starting with melt onset, the melt dynamics 
of the snowpack and the formation and evolution of melt ponds were studied. In addition, the 
physical properties of melt ponds were investigated, mostly in connection to ecological and 
biogeochemical properties. Surface properties and melt ponds were linked to process studies 
on radiative and heat transfer, due to their large impact on the surface energy budget. 
During summer and autumn, the flux of snow and ice meltwater into the ocean played a key 
role with immediate links to the ocean properties, the ecosystem and biogeochemical fluxes. 
The formation and evolution of the meltwater layer in leads and under the sea ice was studied 
in detail in a joint effort across the teams (Section 4.2). 


Dynamics, mechanics, and ridges 

Studies of sea ice dynamics and mechanics were carried out during all legs in the CO. They 
give insights into the connections of sea ice material properties at the millimeter scale to the 
visco-plastic behavior of sea ice at the floe scale. Sea ice deformation at the floe scale was 
measured down to millimeters in the CO using the laser strain instrument. Ridge and lead 
formation were observed through the many dynamical events in the CO by e.g., airborne 
mapping by ALS, the EM-Bird, and aerial imagery. Over 120 position buoys in the DN 
contributed to a regional-scale data set of sea ice drift and displacement. The ice radar system 
onboard of Polarstern also mapped the dynamics within 5 nm around the vessel (Jakel et al., 
2021). 


Ridges and leads had strong logistical implications for the work and fate in/of the ice camp. 
After ridge formation, the temporal development of the consolidated layer, the macro- 
porosity and the keel depth were studied in a number of ridges using the described mass 
balance, ice coring, TLS, transect and ROV methods (Sections 2.3, 2.4, 2.5, 2.10, and 2.11). In- 
situ stresses were measured on the ice floe and will be correlated to atmospheric and oceanic 
forcing as well as the kinematics of ice drift at different scales. 
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The mechanical properties of the ice were studied through in-situ borehole jack testing with 
two different jacks and through sampling ice cores and testing uni-axial strength in the field. 
The full annual cycle of data gives a unique data set that can be correlated with small-scale 
physical properties (temperature, density, salinity) and ice texture allowing for a 
quantification of small-scale (0.01 m) ice mechanical properties and their seasonal variability. 


Microwave interaction with snow and sea ice 

Measurements in the microwave domain (radiometers, radars, reflected GNSS signals) at 
different frequencies and polarizations were conducted simultaneously with the extensive 
physical ice and snow measurement program during the whole MOSAIC drift. The combination 
of all measured frequencies and polarizations allows a better understanding of the processes, 
as well as the volume fraction of emitting and scattering constituents and their geometry, 
which define the seasonal development of microwave signals. Key processes influencing the 
microwave radiation are, e.g., snow metamorphism and ice lenses/layers in the snow, snow 
liquid water content, wicking of brine into snow and desalination of ice during summer, 
changes of ice thickness and the snow/ice temperature profiles (Ulaby and Long, 2014). These 
processes can cause major changes and fluctuations in the microwave signals and can add 
significant uncertainties in satellite sea ice retrieval algorithms. To develop new satellite 
retrievals and quantify uncertainties of existing retrievals, first, the processes that influence 
the microwave retrievals have to be understood better. Second, based on the improved 
process understanding, better microwave emission and scattering models (Tonboe, 2010; 
Picard et al., 2018) can be developed. They require a number of poorly constrained input 
parameters (e.g., correlation length scales or salinity profiles) for which the MOSAIC snow/sea 
ice physics measurements in combination with the microwave measurements are essential. 
The same microwave radiative transfer models (sometimes called forward operators) are also 
needed to assimilate microwave satellite observations in climate and weather prediction 
models. 

The full seasonal cycle was captured by the multi-frequency remote sensing observations. 
Largest uncertainties in satellite retrievals occur when environmental conditions change, like 
during melt-refreeze cycles, warm air intrusions, or rain on snow events, which all were 
observed during MOSAIC and are of special importance. 
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4 First results, case studies and remaining challenges 

This section presents first results from the MOSAIC drift. These results consist of 
methodological advances, which make use of technological developments (Section 4.1) as well 
as scientific results that contrast winter and summer snow and sea ice properties and 
conditions (Section 4.2). The last part (Section 4.3) discusses remaining challenges of the snow 
and sea ice observational program. 


Figure 5 is called the “sea ice clock” and illustrates the sea ice and surface conditions in the 
three COs based on photographs taken by the panorama camera (Section 2.1). Each image 
represents the conditions of one month during the annual cycle: starting from thin and new 
ice in September/October, further ice growth with snow accumulation through the dark 
season until February, into the dense and cold spring ice pack until April, and then snow melt 
and melt pond formation, and the decay of the ice pack in summer. At the beginning of the 
drift, CO1 consisted of SYI with a highly compressed core of the floe and refrozen melt ponds 
over large parts of the floe. Over time new ice formed and the icescape became strongly 
deformed (see below). 


Figure 5: The ‘sea ice clock’. 


4.1 Applying methodological and technological advances 

Methodological and technological advances were implemented for the MOSAIC drift in all 
tasks, allowing a sharp improvement in observing sea ice and its snow cover for an entire year. 
The measurements benefited from development of sensor and instrument technology over 
the last decades, advances in data handling and processing, and improved coordination across 
tasks and disciplines. Technological advances resulted in data sets with higher spatial and 
temporal resolution as well as with reduced uncertainties. It was possible to obtain long time 
series with instrumentation that was earlier only used for case studies or dedicated 
measurements. Covering most of the year, these advanced technologies often turned into a 
backbone of the MOSAIC data set. The technological (including data processing) advances 
enabled the immediate use of measurements to plan additional measurements and aid 
decision making. Selected examples are: 


(1) Combining an eye-safe TLS during winter and a non-eye-safe TLS during summer allowed 
the same regions to be scanned in the CO approximately every 2 weeks (Section 2.5). A custom 
heated enclosure enabled scans throughout the winter in temperatures below —30°C. During 
summer, the Riegl VZ6000 allowed for challenging measurements of wet surfaces. A first 
version of data was processed on board for immediate analysis in support of other 
measurements (e.g., precisely locating the ROV transponders). Figure 6 shows the results of 
two TLS scans of the same piece of ice, containing level ice and a small SYI ridge, on January 
19 and 25, 2020. In this way, the TLS data are used to quantify surface changes related to snow 
fall and re-distribution, as shown in this example. Figure 6d illustrates that the snow 
deposition is concentrated in a snowdrift off the ridge (left of the peak). Similar results on 
snow accumulation in ridged areas are expected from the transect data across the COs, while 
the TLS data provide two-dimensional accumulation details on the centimeter scale. 


Figure 6: Snow redistribution observed via terrestrial laser scanning (TLS). 
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(2) The micro-CT system (Section 2.2) operation in a freezer-lab container allowed for almost 
in-situ imaging of samples, without the delays and likely changes in microstructure. Immediate 
scanning was of particular value for snow samples and the more delicate samples of the 
bottom of sea ice cores, because it ensured that the samples were as close to their original 
state as possible. Sea ice samples from ice cores were transported in insulated customized 
boxes that kept, through eutectic phase change cooling elements, temperatures close to the 
in-situ seawater freezing temperature of approximately -2°C. On the ship the samples were 
centrifuged at this temperature (Weissenberger et al., 1992). This procedure ensured the 
conservation of in-situ microstructure as well as sufficient contrast to retrieve pore network 
details by a subsequent micro-CT imaging (Section 4.3). The immediate onboard visualization 
of the 3D structure of snow and ice gave insights into relevant processes and allowed for 
targeted additional sampling, specifically with respect to the less studied snow-ice interface 
processes. 


(3) The snow pit program was carried out by instrument-based measurements, which resulted 
in a consistent time series with minimal human bias. Main advances resulted from 
standardizing the snow characterization using micro-CT, SMP, NIR and structure from motion 
photography (Section 2.2). These novel quantitative measurement techniques were combined 
with the concept of defining different types of snow pits (e.g., relatively quick snow pits with 
fewer measurements versus more elaborate snow pits using the complete suite of 
measurements; Section 2.2), and thereby ensuring that a standard set of measurements was 
performed for each snow pit type during all legs. 


(4) Studies of dynamical and mechanical properties benefitted from the advance in 
autonomous platforms. In the CO, the laser strain array was operated to investigate sea ice 
deformation and mechanics on the millimeter scale. These data will give new insights into 
mechanical properties, in particular in combination with the measurements of stress and 
strain rates and high precision position monitoring (Section 2.8). 


(5) Optimizing and winterizing the ROV (Figure S3) and surface unit allowed us to operate the 
vehicle in all conditions and without interruptions related to sensors or technology during the 
full annual cycle. Advances in sensor technology allowed for more precise under-ice 
navigation and synchronous recording of all data streams. Specialized additional sensors, e.g., 
the UHI or suction sampling systems, were integrated into the vehicle in the field. In addition, 
the reliable operation of the ROV allowed site selection and instrument maintenance under 
the sea ice. Particularly the sensitive manipulation skills of the ROV system allowed the 
deployment and recovery of sediment traps, inspection and maintenance of under-ice 
sensors, as well as sampling without the need for human divers. Figure 7 shows the light 
transmittance distribution under sea ice in July 2020. Compiling 6 dive days with focus on 
optical measurements makes it possible to quantify the increasing amount of radiation 
penetrating the sea ice as melt progresses, allowing for detailed analysis of the increasing 
spatial variability and the evolution of individual patches with high or low transmittance. The 
overall evolution shows that the initially uni-modal distribution with a modal transmittance of 
0.03 develops into a bi-modal distribution with modes at 0.15 and 0.26. This shift in modes 
represents the evolution from wide-spread wet and melting snow and surface scattering layer 
to a white ice and melt pond pattern. For illustration, see the surface images of June to August 
in Figure 5. 
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Figure 7: Light transmittance through sea ice measured with the ROV. 


(6) The MOSAIC ICE field program strongly benefitted from advances in drone technology. 
Operation of (small) drones for systematic aerial photography and optical measurements (see 
below) advanced the study of surface processes. Advances in and miniaturization of both 
drone and sensors systems allowed for successful operation of these platforms at high 
latitudes, and accurate measurements of broadband and spectral irradiance. The latter was 
achieved in part through the stabilization of sensors using miniaturized gimbal systems that 
could be readily carried by the drones. Additionally, continued advancement of battery 
systems allowed for extended flight times for these platforms, increasing capabilities with 
respect to sampling area and achievable altitude. The drones supported the alignment of 
various measurements at the surface, e.g., the optical measurements along the albedo lines 
(see Figure 9), as well as mapping of the surface at high spatial resolution to provide enhanced 
information over limited spatial coverage within the CO (Section 4.3). 


(7) Data from the ALS and the real-time navigation solution of the inertial navigation system 
of selected flights covering the CO were processed with a short delay of 1 or 2 days after the 
surveys. The processing steps included the estimation of ellipsoidal elevations of the ALS range 
measurements, the transformation of the geographical coordinates into a local cartesian 
coordinate system using position and true heading data from Polarstern as well as the merging 
of all ALS swaths into a single digital elevations model with a spatial resolution of 0.5 m. 
Visualizations of these models served as valuable science planning information in the absence 
of large-scale aerial photography during the winter season (Figures 2 and 8). More than 20 
ALS surveys provided floe maps throughout MOSAIC, with the largest survey consisting of 
three coordinated flights in late December 2019 and covering approximately 70 km2. 


(8) A newly developed sensor suite (Section 2.12) allowed to observe different important 
surface properties at the same time. The topography/roughness of the sea ice cover can be 
directly related to areal properties like ice thickness, melt pond coverage and surface albedo. 
With a lateral resolution of 0.5 m and an elevation uncertainty between 2.5 cm (center) and 
10 cm (edges), the laser scanner data allow digital elevation models to be derived with a 
precision never reached before on airborne operations from Polarstern. The TIR camera has a 
precision of 0.2 K and an accuracy of 1 K and allows detailed mapping of thin ice thickness 
distributions and lead coverage. The digital single lens reflex and hyperspectral cameras 
operated during daylight time allow retrieval of albedo and detailed snow and melt pond 
optical properties. Sea ice thickness surveys (EM-bird) in close temporal vicinity to the 
camera/ALS flights resulted in synergistic measurements that allow a novel areal view of sea 
ice properties. The drones were also used for floe mapping to document and plan the work 
under daylight conditions, thus supplementing the ALS maps. 


(9) The deployment of novel and advanced buoy systems in the CO and DN strengthened the 
acquired time series particularly during the time when Polarstern had to leave the MOSAIC 
floe. In addition to more classical IMBs and Snow Buoys, emerging technologies such as the 
continuous measurements of the vertical profile of sea ice salinity (so called salt harp), in-ice 
and water column light and characterization (optics probe) of under-ice biological activity by 
the means of fluorometry and autonomous acoustic profiling were key contributions of the 
ICE team to the interdisciplinary efforts in the DN. 
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(10) Novel on-ice remote sensing instruments were developed for MOSAIC (Ku/Ka-radar, 
GNSS-R, HUTRAD and L-band microwave radiometers) to evaluate measurements of existing 
and upcoming satellite missions (Section 2.13). For the first time, a large set of 14 remote 
sensing instruments, covering both active and passive measurements as well as a wide range 
of frequencies and polarizations, was operated on the same sea ice floe. Several high- 
resolution optical and synthetic aperture radar (SAR) satellite sensors performed specific 
acquisitions for the MOSAIC region. In particular for SAR, a unique dataset from newly 
launched sensors (e.g., Radarsat Constellation Mission, SAOCOM) and wide range of 
frequencies (ALOS-2, TerraSAR-X, PAZ, COSMO-SkyMed, KOMPSAT-5) were acquired, which 
will allow novel multi-frequency and full polarimetric sea ice analysis at unprecedented 
temporal resolution. These acquisitions required placing an order of the satellite scenes about 
two days in advance on a daily basis. This process was facilitated by the availability of a 
dedicated near-real-time drift forecast product for the MOSAIC floe, to have Polarstern as 
central as possible in the satellite scenes. The product was provided by the Year of Polar 
Prediction (YOPP; Jung et al., 2016) Sea Ice Drift Forecast Experiment (SIDFEx). For MOSAIC a 
consensus ensemble forecast product, based on the different forecast systems, was used. The 
consensus forecasts were provided onboard Polarstern to support decision making. Figure S5 
exemplifies the consensus forecasts product, showing the forecast issued on February 24, 
2020, and revealing that, while the drift over the subsequent 4 months toward Fram Strait 
was on average considerably faster than anticipated, the trajectory remained within the 
ensemble uncertainty margins. 


4.2 Advances through co-location and synchronization 

Coordinating measurements and sampling strategies within the ICE team (across tasks) and 
with other teams allowed the required systematic and representative observations of multiple 
components of the coupled system to be obtained. Here we present examples of such 
activities. 


Probably the most striking and obvious advantage of MOSAiC was the interdisciplinary 
coordination across the scientific teams. Although most field measurements benefitted from 
this advantage, the most prominent example is the coordination of sea ice coring across 
disciplines. Traditionally physical scientists are interested in the ice-snow or ice-atmosphere 
interface and reference their work from the ice surface (i.e., the zero-reference level is set at 
the ice surface). However, most (sea ice) ecological studies are primarily interested in the 
highly porous underside (i.e., bottom) of the ice, where the majority of the biomass is found 
(Manes and Gradinger, 2009). Thus, they tend to reference their work from the ice bottom 
upwards. This difference often makes direct comparisons of data collected with the two 
approaches difficult. During MOSAIC the sectioning of ice cores for physical and biological 
sampling was optimized such that sectioning was initiated both from the top and bottom, and 
the middlemost part of an ice core was left with the “odd” length section. Co-location of coring 
with the ECO and BGC teams captures the seasonal evolution of both physical and biological 
properties for FYI and SYI. Approximately 1500 ice cores were collected during approximately 
30 visits at FYI sites and 25 visits at SYI sites. This coordination effort will allow better direct 
comparison of ice core data across teams. 


Snow sampling was coordinated with other teams, usually collected by the team making snow 
pit measurements, independently of the purpose and later analysis of the snow samples. This 
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strategy ensured that samples from each sampling event were taken at the same site(s) and 
same time and in a consistent and reproducible manner throughout the campaign. Special 
emphasis was given to sampling co-located with other measurements, e.g, at the remote 
sensing site. For interpretation of the remote sensing measurements the snow and ice physics 
measurements are critically needed. 


The observations of low-salinity meltwater layers under the sea ice and in leads during 
summer resulted in a dedicated interdisciplinary program across all teams. Drill lines spanning 
in total 900 m were completed covering different ice types with the specific aim to measure 
the ice thickness and ocean temperature and salinity profiles directly beneath the ice and to 
map the presence of ‘false bottoms’ formed at the meltwater-seawater interface. These 130 
ice thickness observations, of which 30 have associated ocean temperature and salinity 
profiles, overlapped with existing long-term mass balance observations (stakes and thickness 
surveys). Additional melt pond, ice core and under-ice and lead water samples were collected 
by the ICE, ECO and BGC teams to better understand the meltwater sources, and the role of 
these meltwater layers in physical, chemical and biological processes in summer. 


MOSAIC offered the opportunity to co-locate measurements of sea ice surface and bottom 
topography, the internal structure of pressure ridges and the biology associated with ridges. 
This approach allows a complete 3-dimensional view of ridges. Figure 8 shows the surface and 
bottom topography as well as a visual photograph of the pressure ridge called ‘Jaridge’. This 
ridge was located between the SYI closer to Polarstern and the FYI where the ROV and the FYI 
coring site were located (Figure S3). Jaridge was formed during deformation events primarily 
during winter (Leg 2) and included parts of the floe adjacent to the former dark site coring 
site. The Jaridge ridge consisted mostly of 20-40 cm thick ice blocks. This site was used for a 
large suite of coordinated measurements and installations: an IMB was deployed, several 
transects with drilling for ice thickness and consolidated layer thickness were undertaken, sea 
ice cores were collected (ICE and ECO), sediment traps were deployed with the ROV, ROV 
multibeam and the UHI was used to map ridged and level ice, biological samples were taken 
in the ridge, and airborne mapping was conducted with ALS and cameras. This will provide a 
comprehensive and unique view on the evolution of sea ice ridges in the Arctic summer. 


Figure 8: The Jaridge ridge and adjacent lead in summer 2020. 


The evolution of the sea ice surface during the summer is a culmination of changes in the 
radiative energy and mass budget due to variation in snow cover, melt ponds and bare ice. 
Capturing the heterogeneity in surface conditions, in both space and time, is key to 
understanding the seasonal evolution of surface albedo and thus surface energy balance. 
Surface measurements of snow and melt pond properties, and albedo (both spectral and 
broadband), were coordinated with aerial measurements of both surface topography and 
albedo. These observations also overlapped with TLS scans (Section 2.5). Micro-CT 
observations of the ice and snow provide novel insights into the evolution of the surface 
microstructure in relation to the observed radiation budget. In addition to the surface-based 
observations of ice and snow, surface albedo measurements and imagery of the surface were 
collected using drones. These activities included flights over the surface albedo lines to 
compare observations from surface pyranometers and spectral radiometers with those on the 
drones. To align measured albedos along the measurement survey line, drone GPS 
coordinates were corrected to account for ice drift, and the albedo was averaged every 2 m 
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from 4 individual flight passes along the line. Mean drone-based albedo values and standard 
deviation measured along the “lemon drop” albedo line on July 06 are shown in Figure 9 (red 
line). The drone operating along this albedo line carried a modified Kipp and Zonen PR1-V 
broadband pyranometer, measuring wavelengths between 310 and 2700 nm. Surface and 
drone-based measurements show good agreement, with differences primarily resulting from 
the difference in sampling altitude (1 m for surface-based measurements and 15 m for drone- 
based observations), with the higher sensor aggregating the influence of surface type 
variability. This influence of sensor altitude on observed surface albedo is being investigated 
in a separate study. 


Figure 9: Surface albedo during summer. 


Figure 10 shows winter and summer time-series measurements of several co-located remote 
sensing instruments (Table 2 and Figure S4). During a storm event in November 2019 
measurements by both radars and radiometers show changes in backscatter and brightness 
temperature, respectively. Due to the fact that air and surface temperatures reach a 
maximum of —5°C, melting plays no role in these changes, but snow metamorphism already 
can happen. Also, strong winds of up to 20 m s* caused snow compaction and redistribution. 
The radars and radiometers show stronger changes for higher frequencies and partly opposite 
behavior depending on frequency. The 35 GHz radar backscatter decreases by close to 2 dB 
while at 15 GHz it increases by about 1 dB during the event. At the lowest frequency of 1.3 GHz 
no changed radar response is observed during the event, which is in correspondence with the 
minimal changes in GNSS-R reflectivity operating at a similar frequency. During the storm 
event the increase in magnitude of microwave brightness temperatures (TB) is increasing with 
frequency: e.g., at 89 GHz TB is increasing by about 50K, which corresponds to a strong 
emissivity change, while at 1.4 GHz almost no TB increase is observed. The peak of the storm 
on November 16, 2019, could not be captured by most of the remote sensing instruments 
because a crack opening, due to strong ice dynamics, at the site caused a power outage. 

In summer, on September 13, 2020, even stronger fluctuations in both backscatter and 
brightness temperatures occurred when the air temperature became positive and melting 
increased. Rain on snow additionally changed surface properties. GNSS-R reflectivity, both at 
the remote sensing site and on Polarstern, increased during the event. On September 15, after 
the temperature dropped below zero again, masurements also stabilized. Similar to the winter 
event, the changes in magnitude and variability increased for higher frequencies both for the 
radars as well as the radiometers, with the exception of the 1.4 GHz radiometer. However, 
changes were much stronger in summer compared to winter (note the different y-scales in 
Figure 10). 

Ice thickness influences the signal at very low microwave frequencies. Measurements by the 
1.4 GHz L-band show lower brightness temperatures at a partly refrozen lead (light blue line) 
than for the thicker ice at the remote sensing site (dark blue line; Figure 10). This effect is used 
to retrieve thin ice thickness from L-band satellite measurements from SMOS and SMAP 
(Kaleschke et al., 2016; Patilea et al., 2019) and in future CIMR. The higher measurement 
variability at the lead likely is due to the changing ice conditions caused by ice dynamics. 
Variability in brightness temperature and backscatter induced by environmental changes like 
those shown in these two case studies are not accounted for in current satellite sea ice 
retrievals for ice concentration, thickness, or snow depth and are causing uncertainties in the 
satellite data sets. In conclusion, during both events lower frequencies were less affected by 
the environmental changes and thus should result in more stable satellite retrievals. By 
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combining the on-ice remote sensing measurements with the detailed snow and ice physics 
measurements, we can better understand and model the interactions of microwaves with 
snow and ice. Future work will improve satellite remote sensing methods to obtain better sea 
ice satellite climate records. 


Figure 10: Remote sensing signatures during winter and summer. 


4.3 Case studies in winter and summer 

To demonstrate the comprehensiveness of the snow and sea ice measurements during 
MOSAIC, two weeks of measurements are selected as case studies. The week starting January 
20, 2020, (Leg 2) represents the work during winter darkness in the polar night and the week 
starting on July 06, 2020, (Leg 4) represents the work during summer. Figure S1 shows how 
the work from the different tasks was distributed over the respective week. Additional optical 
measurements were carried out under daylight conditions. 


Sea ice, snow and weather conditions differed significantly for both case study weeks. As 
expected for high Arctic conditions, the mean air temperature was —29.0°C during the week 
in January and 0.3°C during the week in July. Wind speed at Met City (Figure 2) was 4.6 ms? 
and 6.4 ms for the weeks, respectively. In general, the winter period was characterized by a 
transition from weak winds with periodic thin, ice-dominated clouds to a couple of passing 
snow storms later in the week with increased winds. The summer week was similar in some 
regards, but with much higher temperatures, characterized by persistent low-level, liquid- 
dominated cloud and fog, with two short storms that produced light drizzle near the surface. 
A detailed description of the meteorological conditions at the CO is in preparation led by M. 
Shupe. 


To contrast the different snow and sea ice conditions, Figure 11 illustrates various aspects of 
the surface and ice conditions. The figure shows the CO and its surrounding on composites 
from hundreds of aerial photographs taken from the Polarstern helicopters. During polar 
night, TIR photos were used to map the surface temperature, which indicates sea ice types 
based on differences in surface temperature. The map shows a dense ice pack with few 
deformation zones and elongated structures of thinner ice (higher temperatures) resulting 
from leads with new thin sea ice. During summer, visible photos were taken and give a view 
of the surface conditions. Towards the end of June and in July, the sea ice surface consisted 
of many distinct floes, widely covered with melt ponds. The observation over the year of the 
same ice pack and area will allow consideration of how much of the wintertime surface and 
structural features determine summer features, such as the melt pond distribution. In general, 
the observational concept of TIR and ALS flights during the dark season proved most useful 
(Figure 2). Both methods were used to map surface conditions and guide the on-ice 
measurement program as well as to support logistics through results on snow and sea ice 
surface properties, immediately processed on board Polarstern. 


Figure 11: Aerial photo mosaics of the central observatory and its surroundings. 
The drift and deformation of the ice pack was one of the most remarkable features of the 
MOSAIC field experiment, and more extensive than expected. In total, the ice camps covered 


distances of 2,354 km (Drift1 with CO1), 622 km (Drift2 with CO2) and 345 km (Drift3 with 
CO3) based on hourly position data during the manned observations. These tracks are shown 
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in Figure 1, which also shows the continued drift of the COs after Polarstern has departed 
(Table 1). Satellite data indicate that the MOSAIC drift was about 25% faster than the 
climatological mean drift, and was caused by large-scale low-pressure anomalies that 
prevailed around the Barents-Kara-Laptev Sea region between January and March (Krumpen 
et al., 2021). Ice drift and deformation was tracked in the DN as well as in the CO. Figure 12 
compares the ice drift speed with the wind speed for January and July and discriminates 
whether the floe was in free drift or moved with the closed ice pack. It reveals much higher 
drift speeds with a mean of 0.19 m sin July, when the floe was in free drift, while drift speeds 
with a mean of 0.10 m s‘twere only about half that value in January, when the floe was 
embedded in the interior ice pack with other floes. In addition, the variability of drift speed 
was 1.5 times higher during free drift in July than in January. Additional analyses of this data 
set for the entire year will reveal details of the seasonality of different drift modes as well as 
dominating processes, as the start of inertial motions became more important around mid- 
June (Dethloff et al., 2021). Seasonal contrasts in drift and deformation patterns can also be 
seen in the ice radar animations over the entire CO1 (see reference to video animation in data 
section) Future work will analyze the interaction between ice fractures and wind forcing also 
on larger scales by combining the wind stress information across the DN. 


Figure 12: Sea ice drift speed as a function of mean wind speed. 


Figure 13 shows the rate of deformation (,/divergence? + shear”) derived from two 


consecutive Sentinel-1 SAR images in a 200 x 200 km distance centered around the MOSAIC 
CO for winter (December 2019) and summer (June 2020) conditions (for method see Albedyll 
et al., 2021). Intersecting lines of strong deformation (linear kinematic features) are present 
both in winter and summer ice pack. Differences in winter and summer deformation, e.g., in 
the degree of localization of the deformation visualized in the width of the linear kinematic 
features (Figure 13), are caused by changes in the ice pack strength. The seasonally varying 
deformation provides additional insights to the different sea ice conditions previously 
presented in Figure 11, and adds the larger scale to the airborne observational data. 
Additional analyses of spaceborne remote sensing are described in Krumpen et al. (2021). The 
SAR-based ice deformation can, in future, be combined with the buoy observations from the 
DN and the ship-radar-based deformation to obtain data products of high spatial and temporal 
resolution. Thus, these observations may be used towards an improved understanding of ice 
dynamics on different spatial scales. 


Figure 13: Total deformation derived from two consecutive Sentinel-1 SAR images. 


On the floe scale, Figure 14 describes the distribution of sea ice and snow thickness for January 
and July, and melt pond depth (for July only). The transect data (Section 2.10) shows a winter 
ice mode of 1.3 m and a thicker summer mode of 1.8 m (mean 2.7 m). The winter and summer 
modes of snow and sea ice thickness distributions demonstrate the effects of thermodynamic 
ice growth/decay and deformation throughout the seasonal cycle. The winter distribution 
features two distinct total (snow plus ice) thickness modes associated with a thinner young 
ice and a thicker mix of level/deformed FYI and SYI. In summer, in turn, the distribution 
exhibits two modes of level FYI and SYI accompanied by a broad tail of thicker ice formed by 
deformation processes during the preceding winter. The transect lines between winter and 
summer had only partial overlap. As one would expect, snow thickness was much reduced in 
summer compared to winter/spring, which is obvious from the reduction in modal snow 
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thickness from 0.06 to 0.03 m. In July the modal melt pond depth was 0.1 m, with maximum 
depth of 0.5 m. The large coverage of melt ponds in summer is shown in Figure 11b. Results 
from a helicopter electromagnetic ice thickness survey on July 01, 2020, (data not shown) 
revealed a thicker summer modal total thickness of 2.1 m, but a thinner mean total sea ice 
thickness of 2.5 m over the CO and the surrounding ice within a radius of 50 km. 


Figure 14: Sea ice thickness and snow depth distributions. 


The characteristics and evolution of FYI is exemplified in Figure 15 and data are given in Tables 
S2 and S3. In winter (December) the upper 15 cm from the ice surface consisted of granular 
ice (Figure 15a), while most of the core (19-65 cm) was composed of columnar ice. A 6-cm 
thick skeletal layer was apparent at the core bottom (65—71 cm) of the growing FYI. While the 
ice in December (Figure 15b) had the typical C-shape salinity profile and was cold due to the 
low atmospheric temperatures, the rapid warming and desalination of FYI in summer is shown 
with the data from early July (Figure 15c). 


Figure 15: Sea ice physical properties. 


The surface properties of the topmost layer of sea ice showed large temporal and spatial 
variations throughout the MOSAIC observation period. The typical winter snowpack consisted 
of an often wind-packed new snow layer at the surface, with low density, large specific surface 
area and relatively small optical equivalent diameter (0.00—-0.15 mm; Figure 16a). Due to the 
large temperature difference between the air and the ocean, the snowpack experienced large 
temperature gradients during most of the winter, which caused strong recrystallization and 
faceting into columnar snow structures with a strong geometric anisotropy (depth hoar). We 
observed the formation of euhedral crystals in the low-density regions (0.40—0.65 mm; Figure 
16a), in contrast to subhedral crystals in the layers with higher densities (0.25-0.40 mm; 
Figure 16a). Our data show that the strong temperature gradients in the Arctic snowpack 
during winter cause a highly anisotropic microstructure that has significant effects on the 
thermal conductivity of snow and sea ice growth. In addition, they affect microwave 
properties for remote sensing (Figure 10). The sea ice surface in summer (Figure 16b) is 
generally free from snow and consists of the SSL that is formed by melting and draining 
processes in the sea ice. In the field visually distinguishing between snow and the SSL is hard; 
however, data from the micro-CT highlight the microstructural differences between the SSL 
and snow (Figure 16b). The density and the optical equivalent diameter increased with depth, 
reflecting the transition from the drained and melted sea ice at the surface to values typical 
of frozen sea ice (approximately 700 kg m°3) toward the bottom of our sample. The specific 
surface area was relatively consistent, showing an increase at the surface due to preferential 
melting of columnar ice. The strong anisotropy of the SSL reflects the pre-existing structure of 
the original columnar sea ice crystals and brine channels and influences optical properties of 
the sea ice surface as well as remote sensing retrievals. 


Figure 16: Exemplary vertical profiles of the sea ice surface and snow during winter and 
summer. 
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4.4 Remaining challenges of the snow and sea ice program 

The observational program of MOSAIC had to compromise in various aspects. The largest 
impact on the observations of the full annual cycle was that Polarstern had to leave the main 
floe (CO1) temporarily to exchange personnel and to supply the vessel. This absence resulted 
in a gap in manned observations between May 16 and June 19, 2020, (Table 1), a key phase 
of the annual cycle: the spring-summer transition and the early melt onset. The formation 
and decay of sea ice during the shoulder seasons is of increasing importance in the increasingly 
seasonal ice cover of the Arctic. Studying these processes will need dedicated expeditions to 
observe sea ice and snow processes of thin and very fragile ice covers in very close connection 
to oceanographic and atmospheric conditions. During MOSAIC, only autonomous 
measurements in the DN were able to bridge this gap to some extent and to obtain some key 
parameters continuously across this gap. Nevertheless, further studies of similar 
interdisciplinary complexity are needed to add more coordinated in-situ observations of this 
phase, where small changes in timing have large impacts on the total energy budget and the 
seasonality of ecological processes (Nicolaus et al., 2012). 


All observations are limited to one specific drift trajectory along the Transpolar Drift from 2019 
to 2020. This leaves the challenge of upscaling and generalization, which will mostly rely on 
merging the in-situ observations with numerical models (Section 5.2) and remote sensing 
(Section 5.3). Although many relevant sea ice processes were very well covered, individual 
parts are likely missing or lacking details. The work on the snow and sea ice properties in the 
DN was, in the end, quite limited both in space and time. The need for a permanent ice camp 
and fully manned research icebreaker for most aspects of the field program excluded most 
observations of marginal ice zone processes. 
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5 Linkages within the coupled Arctic system 

MOSAIC demonstrated the benefit of improved coordination across disciplines, which allowed 
us to study linkages and to realize ambitious plans in integrating methods and disciplines. 
Section 5.1 describes linkages to the individual sub-systems (including the work of the other 
four teams) and points to the comprehensive data sets, which will build the legacy of MOSAIC. 
The sampling strategy was guided by the needs of sea ice and climate models as well as for 
satellite remote sensing. These linkages are discussed in Sections 5.2 and 5.3, also in 
connection to the DN, because they allow upscaling of observations relevant for climate 
models and remote sensing. 


5.1 Linking the individual sub-systems 

Combined analyses of key parameters across all sub-systems will allow unprecedented studies 
of the coupled Arctic system. The physical properties of snow and sea ice are crucial for 
interpreting biogeochemical and biological observations as well as for understanding 
atmospheric composition and aerosols (Figure 4). For example, the evolution of snow and sea 
ice microstructure (Figures 14 and 15), porosity and permeability, is known to affect the 
transfer of trace gases (e.g., Nomura et al., 2018) and aerosol particles across the atmosphere- 
ice-ocean boundary (Frey et al., 2020). The snow on sea ice and the sea ice both act as a 
chemical reservoir and reactor releasing a range of chemical trace gases, aerosol particles and 
their precursors with significant impacts on tropospheric composition and oxidizing capacity, 
and potentially on climate via influencing clouds (Grannas et al., 2007). Observations of gas 
and particle fluxes across the atmosphere-snow-ice-ocean boundary were carried out at 
different scales by the BGC team and ATMOS team. Here both snow and sea ice properties 
from different ice types as well as ice dynamics (especially formation of leads, e.g., Figure 13) 
are likely crucial factors for the observed fluxes. The close coupling of the snow and sea ice 
cover with atmospheric forcing (especially wind forcing) became obvious (Graham et al., 
2017). The dynamics of the ice pack challenged the field program, but dynamics are critical for 
a better understanding of energy and mass exchanges between atmosphere, sea ice, and 
ocean. Here we gave examples of dynamical processes with respect to the fast drift (Figure 
12), the formation of leads and ridges (Figures 8 and 13), and their effects on the sea ice 
thickness distribution (Figure 14). While the sampling activity (coring, snow sampling) 
captures the seasonality, more rapid or transient changes (e.g., due to warm air advection) 
are more reliably covered by continuous measurements of autonomous systems; their 
combination provides new insights into the coupling between snow and sea ice physical 
properties and atmosphere—ocean exchanges. 


The role of melt ponds and open water fraction was investigated intensively with respect to 
formation processes (Figure 5), energy budgets (Figures 7 and 9) and the impacts of their 
large-scale distribution (Figure 11) on the sea ice mass balance (Figure 14). Findings from these 
investigations will also link to concurrent ecological and biogeochemical studies. The small- 
scale surface heterogeneity and, in particular, the contrasts of surface temperature, trigger 
exchange processes between the Arctic system components. During low ice surface 
temperatures, the warmer open water areas or partly refrozen leads act as a source of latent 
and sensible heat, which modifies the atmospheric boundary layer and triggers atmospheric 
convection (e.g., Lupkes et al., 2008; Schmale et al., 2021). The strong release of water vapor 
into the atmosphere in summer results in enhanced formation of low-level clouds and fog 
(Tjernstrom et al., 2012). While the leads represent a conduit for strong heat loss (Figure 11) 
and new ice formation during winter (Itkin et al., 2018), they are windows in the sunlit 
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summer (Figure 7) period that allow more solar heating of the upper ocean (Taskjelle et al., 
2017) and intensified pelagic primary production before melt ponds appear (Assmy et al., 
2017). 


The observations of sea ice mass balance provide an integrated measure of the energy 
balance, which is a result of the interactions and energy exchange with the atmosphere and 
the ocean. Here both thermodynamic and dynamic forcing are relevant. The first observations 
of platelet ice in the central Arctic during winter (Katlein et al., 2020) suggest that this 
mechanism can also contribute to the mass balance of Arctic sea ice. The detailed 
observations of ocean heat and freshwater content and fluxes by the OCEAN team, provide 
data sets to examine the interaction between upper ocean properties and sea ice under 
different ice regimes. 


The detailed work on ridges and leads (Figure 8), in close coordination and collaboration with 
oceanographic and ecological aspects, represents a leap forward from simply investigating the 
level and modal snow and sea ice properties towards a much more comprehensive 
understanding of the whole ice pack. These measurements allow specific analyses of 
deformed and ridged ice (Gradinger et al., 2010; Fernandez-Méndez et al., 2018; Shestov et 
al., 2018), e.g., through comparisons of sea ice habitats, energy and mass fluxes and 
contrasting processes of level ice with deformed ice, thin ice, and leads. 


5.2 Relation to modelling 
A synthesis of sea ice field measurements with results from numerical models involved in 
MOSAIC is expected to demonstrate how process models with different levels of complexity 
and coupling can help to advance our understanding and prediction of Arctic climate change 
and contribute to improved numerical weather and sea ice prediction. Some of the expected 
progress is related to: 
e the assessment of unresolved, missing or overly simplified or under-represented 
processes in the sea ice component of the fully coupled Arctic system; 
e the exploration of inadequacies and possible improvements along the model 
atmosphere-sea ice and ocean-sea ice coupling channels; 
e the analysis of model sensitivities to the scale-aware sea ice parameter space; 
e the placement of the MOSAIC sea ice measurements within a broader spatial (i.e., pan- 
Arctic) and temporal (i.e., interannual to multi-decadal) context. 


Major challenges remain in how to optimally utilize observational data in different models, 
including those related to 

e using Lagrangian point-observations for evaluation of gridded, space- and time- 
discrete model output; 

e post-processing heterogeneous, erroneous and inconsistent data sets into data 
products that may be useful for model initialization, forcing and evaluation; 

e bridging different scales of in-situ observations (e.g., from the CO) and remotely 
sensed measurements in order to up- or down-scale such data to a model grid cell or 
mesh; and 

e homogenizing the heterogenous data sets of individual parameters into a common 
data product, e.g., by interpolating to common time stamps or by classifying data 
according to the snow and sea ice states. 
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On the model side, future work will be needed to: 

e evaluate sea ice models and their existing parameterizations in a process-oriented 
manner that benefits from an improved process understanding; 

e use process models, large eddy simulations or direct numerical simulations to derive, 
test and make available advanced parameterizations of processes impacting snow and 
sea ice states and variability; 

e assimilate sea ice data to produce a better state estimation and forecast initial 
conditions for near-term forecasts; and 

e derive improved sea ice forecasts and projections across timescales, including the 
consideration of internal and forced variability. 

During the field phase of MOSAIC, different observations were realized to enable a process- 
oriented model evaluation and thus to assess and improve systematic model biases through 
improved process understanding and related parameterizations (e.g., Figure S5). 


From the sea ice perspective, the most immediate connections, and thus model developments 
are currently expected from: 

(1) Improvements to the representation of snow on sea ice (Figures 6, 14 and 16). The 
evolution of the snow cover (fraction, depth, distribution, thermal properties) is often overly 
simplistic in various types of models, e.g., in global coupled climate models (Chen et al., 2021; 
Webster et al., 2021) and models to reproduce re-analyses (Sato and Inoue, 2018; Batrak and 
Muller, 2019; Wang et al., 2019). Model sensitivity studies also indicate that the simulated 
climate is very sensitive to the snow representation (Urrego-Blanco et al., 2016). The new field 
observations cover all relevant processes (Figure 4) from freeze-up to melt, including remote 
sensing observations on large scales. Particular improvements are expected for estimates of 
snow accumulation and re-distribution over deformed ice (Liston et al., 2018), surface 
topography (Figures 6 and 8), thermal properties and heat transfer between atmosphere and 
ocean, as well as internal snow processes (Wever et al., 2020). 

(2) Improvements in surface albedo and optical properties (Figures 7 and 9). The evolution of 
optical properties, along with other energy transfer terms, has huge implications for the total 
energy budget and is still one of the most critical aspects of sea ice modelling and tuning (e.g., 
(Holland et al., 2012; Jakel et al., 2019). Model studies indicate that the treatment of the 
albedo, including factors like melt ponds, modify the transient climate evolution (e.g., Holland 
et al., 2012). They also show the need for inclusion of specific surface albedo for different 
surface types and the factors that drive the evolution of those surface types over time. Based 
on the MOSAIC observational data sets, better implementation of albedo schemes and melt 
pond thermodynamics is expected, leading to improved heat fluxes over the heterogeneous 
pack ice and more realistic surface energy budgets. These improvements will also have 
important implications for atmospheric processes and how they are represented in models, 
due to the strong coupling between surface fluxes and boundary layer processes. 

(3) Improvements in sea ice dynamics (Figures 12 and 13). An improved inclusion of dynamical 
processes, including deformation and sea ice surface/bottom roughness (Figure 8) is critical 
for an accurate description of momentum fluxes across the atmosphere-ice-ocean interfaces 
and internal stress (Hutchings et al., 2011). Here, the formation and decay of ice pressure 
ridges and the sea ice floe size distribution are of particular importance to describe the 
atmosphere-ice and ice-ocean drag and internal ice stress. Both thermodynamical and 
dynamical processes are important for sea ice change (Dethloff et al., 2021). The main benefit 
of the MOSAIC observations is expected from the dedicated program on pressure ridges and 
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leads in the CO (local scale) and the unique DN, which allows deriving dynamic parameters on 
a broad range of scales (e.g., ice radar, helicopter surveys, satellite remote sensing). 


To enable model evaluation and development of improved parameterizations, a main goal is 
to merge the existing observations into an explicit “MOSAiC standard forcing and 
benchmarking data set” for thermodynamic/dynamic sea ice models. A second goal is to 
advance coupling with ecological and biogeochemical processes, where the co-located 
measurements will be most important. They bridge across geophysical and biological aspects 
of snow and sea ice and allow for benchmarking simulations of causal relationships between 
environmental conditions and responses of different components of the Arctic sea ice 
ecosystem. 


Over the last decades, assimilation of sea ice and snow data became an important tool to 
study sea ice and its role in the coupled climate system and to initialize forecasts. The field 
data will be used to improve the assimilation of sea ice concentration, sea ice thickness and 
snow depth in coupled ocean-sea ice models (Fritzner et al., 2019), in particular by calibrating 
and validating satellite-derived products and thus reducing uncertainties in the large-scale 
assimilation data sets (Figure 10). The impact of assimilating sub-grid scale sea ice thickness 
distribution derived from Cryosat-2 and ICESat-2 observations will be assessed and support 
simulations with a stand-alone sea ice model. 


The field phase of MOSAIC also demonstrated how model applications can support and guide 
the field experiment. Sea ice simulations and in particular ice drift forecasts were provided to 
Polarstern in near real time through the SIDFEx initiative (Figure $5). Supplementing the 
atmospheric forecasts, this supported the on-ice measurement program, e.g., in planning 
intensified observation periods. These data sets were mainly used for near-real-time 
verification of the model performance and for real time educational applications. 


5.3 Relation to remote sensing by aircraft and satellites 

As part of the MOSAIC expedition, several airborne surveys were conducted with the AWI 
research aircraft Polar-5 and Polar-6 from Longyearbyen (Spitsbergen) between August 17 and 
September 17, 2020. The aircrafts are of type Basler BT-67 and were operated by the Canadian 
Company Kenn Borek Air Ltd. Calgary. The main goal of these survey flights was to extend the 
Spatially limited observations of the atmosphere and sea ice in the vicinity of CO to a larger 
area. While the airborne remote sensing from the Polarstern helicopter cover scales up to 
100 km with a high, weekly temporal sampling, the aircrafts cover several hundreds of 
kilometers and can carry more instrumentation. When the airborne surveys began in August 
2020, Polarstern was already on its way to the central Arctic to relocate to CO3 (Figure 1), and 
thus out of aircraft range. However, various DN buoys, positioned between 79°N and 81°N in 
early September, were still active and within range. These buoys allowed the sea ice 
conditions over the DN area to be recorded after Polarstern had departed. The surveys also 
complement earlier sea ice surveys made in Fram Strait between 2001 and 2018 (Belter et al., 
2021). For this purpose, the aircraft Polar 6 was equipped with a number of sensors, including 
the EM-Bird, a laser scanner and optical instruments. Together with the laser scanner, the 
optical camera recorded melt pond distribution, surface elevation, floe size distribution, and 
other surface properties. An insight into the different sensor systems and their specifications 
is given in Herber et al. (2021). The track of a flight over the DN made on September 02, 2020, 
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is shown in Figure 1. Average ice thickness in the vicinity of the DN was 1.44 m, while the 
modal ice thickness was 0.93 m. Note that in parallel to the MOSAIC airborne campaign, sea 
ice and oceanographic surveys were carried out by the Norwegian research vessel Kronprins 
Haakon in Fram Strait, and a sea ice floe next to the vessel was overflown while in-situ 
observations were taken. The second aircraft Polar 5 was equipped with in-situ, remote 
sensing, and basic meteorological instrumentation. The observations focused on 
characterizing the Arctic atmospheric boundary layer in conjunction with ocean-atmosphere 
interactions, clouds, solar and terrestrial radiation, and aerosols. More details about the 
airborne atmospheric survey program may be found in Herber et al. (2021). 


Satellite observations can extend the local MOSAIC observations to regional and hemispheric 
scales. Long-term satellite time series and climate data records can put the MOSAIC 
observations into a temporal context and answer questions about how representative or 
unusual was the MOSAIC year (Dethloff et al., 2021; Krumpen et al., 2021). Routinely taken 
satellite acquisitions already provide a large collection of daily sea ice related properties from 
space, which will help and extend the MOSAIC data analysis: sea ice area, thickness, ice type, 
drift, albedo, snow depth, melt pond coverage and more. However, in addition a large 
collection of satellite data was acquired specifically for MOSAIC, especially high-resolution SAR 
and optical data. During sunlight, observations of various optical sensors with different 
footprints are available. WorldView-2 (meter footprint), Sentinel-2 (tens of meter footprint), 
Sentinel-3 and MODIS (hundreds of meter footprints) allow the scaling of in-situ to satellite 
observations at these scales. Several space agencies helped to acquire a unique dataset of 
radar backscatter images at different spatial resolution, frequencies, and polarizations from 
several different SAR sensors (ALOS-2, COSMO-SkyMed, KOMPSAT-5, PAZ, Radarsat 
Constellation Mission, SAOCOM, Sentinel-1, TerraSAR-X). This help allowed monitoring the CO 
and DN and an unprecedented temporal resolution of several acquisitions per day. However, 
MOSAIC was drifting through the “pole hole” of many satellite sensors (Figure 1) for a 
significant amount of time and thus not all observations are available year-around. Retrieved 
quantities from these higher resolution satellites acquisitions will be, e.g., floe size 
distributions, lead and ridge locations, ice types, ice deformation, melt pond and albedo 
distributions. On the other hand, the MOSAIC in-situ observations will help to improve and 
develop new satellite retrievals (descriptions in Sections 2 and 3 and examples in Section 4). 


6 Conclusions 

The snow and sea ice field program was designed to monitor all key parameters of the snow 
and ice system over one full annual cycle, while the intensity of observations over the year 
were adapted with respect to (1) deciphering sea ice evolution from autumn freeze-up to 
summer melt, (2) capturing key events and seasonal changes, and (3) accommodating 
targeted research on emerging science questions. The resulting data set will allow to better 
quantify the causes and consequences of the evolving and diminishing Arctic sea ice cover, as 
one of the central elements of the Arctic climate system. The first results show the active 
dynamics of the entire Arctic ice pack. MOSAiC experienced and observed rapid ice 
transformations and motions, and strong deformations along the drift in all seasons. These 
observations demonstrate the importance of dynamic processes in comparison to 
thermodynamic processes, as well as the role of ridges (deformed sea ice) and leads in 
comparison to level sea ice. Leads, in particular, link the ice pack to atmosphere-ocean 
exchange. The program also captured the evolution of the snow cover, across different ice 
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types, including leads and deformed (ridged) ice. Strong temperature gradients during winter 
cause a highly anisotropic microstructure with significant effects on the thermal conductivity 
of snow. The large spatial and temporal variabilities of the snow pack impact sea ice growth 
and have to be considered in more detail, both in observations and in simulations. These 
details will allow better understanding of feedback processes between the ice and the 
atmosphere, including gas and particle exchange. The snow and ice measurements combined 
with the remote sensing observations will lead to a better understanding of sea ice microwave 
emission and scattering, including their temporal variability which is needed to improve 
satellite datasets. In particular, the role of events like the passage of low-pressure systems 
(storms), warm air intrusions, rain and snow fall, lead and ridge formations can be studied 
using integrated data sets from all key parameters. Ice-ocean interface processes, such as the 
formation of platelet ice, the summer freshwater layer or the roughness of the ice-ocean 
interface, were observed in great detail and will likely guide upcoming research with respect 
to the changing Arctic sea ice cover. 


Even though the MOSAIC snow and sea ice work was extremely comprehensive, some aspects 
could only be studied in a very limited way. For example, MOSAIC was only marginally able to 
study the transitions of sea ice from/into open water. Dedicated studies in the marginal ice 
zone may be able to shed more light on this topic in the future. Furthermore, the temporary 
departure of Polarstern from the ice due to logistical reasons resulted in a gap in many of the 
in-situ observations during melt onset and transition into summer. Despite these limitations, 
we expect that the overall work will lead to a better process understanding of snow and sea 
ice and their linkages to atmosphere and ocean as well as improved forecast capabilities. Our 
expectation is that this work will strengthen research on the global coupled climate system, 
in particular with respect to the seasonally frozen Arctic Ocean. 
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Figure legends 


Figure 1. Drift tracks of the central observatories (CO) of MOSAIC in 2019-2020. 

Colored segments illustrate the month of the drift of the COs (Table 1) with thick parts 
indicating the manned drifts, and thin parts the unmanned drift. The black lines show the 
tracks of Polarstern. The grey line shows the drift track of the floe prior to the MOSAIC drift 
and the dashed black line, the flight track of the MOSAIC aircraft campaign on September 02, 
2020. Circles indicate the pole holes of different satellites: 89°N, e.g., AMSR2; 88°N, e.g., 
CryoSat2, IceSat-2, SMOS; 87.5°N, e.g., Sentinel-1. The background shows the sea ice extent 
for the March 2020 maximum and the September minima for 2019 and 2020. Labels denote 
the start of the legs and COs. 


Figure 2. Main sites and installations in the central observatory at the end of Leg 3 (CO1). 
Note that panel B extends the map in A to the left (join the arrows). The background shows 
the airborne laser scanner image from April 23, 2020 with grey areas indicating “no data”. Due 
to active ice dynamics, the positions of some sites were approximated, and some sites were 
outside of the map range or were destroyed. The term “old” refers to sites that were active 
earlier during the expedition but were no longer maintained at the end of Leg 3. 


Figure 3. Snow and sea ice observations during the field phase of MOSAIC. 

Each row shows dates when continuous (lines) or discrete (points) measurements were 
performed. Dark grey and white areas indicate the legs, light blue areas indicate the time 
when Polarstern drifted with the COs, light grey columns (during Legs 2 and 4) indicate the 
times of the case studies described in Section 4.3. (A) On-ice and airborne measurements, (B) 
remote sensing instruments on the ice and on Polarstern. Abbreviations: Hyperspectral 
camera, IR/video camera. Acronyms are listed in Table $1. Color description for rows with 
multiple colors from top to bottom (red: r, green: g, blue: b, magenta: m, yellow: y): 

Drone flights: (r) HELIX, (g) Mavic and Spectra Quadrocopter; 

Ice mechanics: (r) Hydrocomplex, (g) bore hole jack, (b) both; 

Albedo lines: (r) spectral + broadband, (g) spectral only, (b) broadband only; 

Ponds: (r) remotely controlled water color spectroradiometers, (g) hand-held water color 
spectroradiometers, (b) water sample, (m) pond depth; 

Coring: (r) FYI, (g) SYI, (b) FYI+SYI, (m) Ridge, (y) other sites; 

ROV: (r) Sensor data only, (g) Additional net on ROV; 

Helicopter: (r) Laserscanner + cameras, (g) HEM bird, (b) Helipod; 

Bridge observations: (r) at station, (b) during transit; 

Panomax: (r) daylight conditions, (b) darkness; 

Radiation station: (r) daylight conditions, (b) darkness; 

C-SCAT: (r) regular measurements, (b) recorded data with problems; 

HUTRAD: (r) regular measurements, (b) recorded data with problems. 

ARIEL: (r) regular measurements, (b) recorded data with problems. 


Figure 4. Schematics of the snow and sea ice processes studied during MOSAIC. 
The close-up illustrations focus on the snow-ice (left) and ice-ocean (right) interfaces. 


Figure 5: The ‘sea ice clock’. 
Photographs from the MOSAIC floe arranged clockwise with one picture per month (January 
to December). Photographs were recorded by the panorama camera on the crow’s nest of 
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Polarstern and show the view ahead of the vessel. Photos were selected to be representative 
for the month while providing best available image quality. 


Figure 6. Snow redistribution observed via terrestrial laser scanning (TLS). 

(A) Complete TLS scan from January 25, 2020, colored by topographic relief, (B) 40 x 40 m 
close up of that scan, (C) perspective view from the viewpoint (eyeball and yellow triangle in 
(B) towards the ridge, and (D) comparison of the topographic relief along the transect line for 
both scans, with 95% confidence intervals for the surface reconstruction in each profile shown 
in shading. 


Figure 7. Light transmittance through sea ice measured with the remotely operated vehicle. 
(Top row) transmittance (integrated 350-920 nm) along the dive track at 2 m depth for 6 dives 
in July 2020. (Bottom row) histogram of light transmittance derived from the above dive. 


Figure 8. The Jaridge ridge and adjacent lead in summer 2020. 

(A) Approximate surface elevation of the sea ice or snow surface from airborne laser scanning 
on June 30, 2020, (B) Bottom topography from the multi-beam sonar on the remotely 
operated vehicle (ROV) Beast on June 25, 2020, (C) Surface photography on July 04, 2020, and 
(D) Surface photography with surface installations and Polarstern in the background on July 
21, 2020. In B) the location of the ROV hut, sediment trap deployment, IMB (2020M26) 
deployment and the approximate views of C and D are indicated. The dashed lines indicate 
the locations of consecutive drilling transects across the ridge. 


Figure 9. Surface albedo during summer (Leg 4). 

(A) Broadband surface albedo and (B) spectral surface albedo along the 200-m-long Lemon 
Drop line on July 06, 2020, (Leg 4). Broadband measurements in blue and spectral 
measurements were completed on the ground (1 m height), and broadband measurements 
in red were taken from the drone (HELIiX) with a flight altitude of 15 m. Shading represents 
the standard deviation of 2-m spatially averaged measurements (recorded at 1 Hz). (C) 
Photomosaic taken from the drone with the red line showing location of the drone measured 
albedo line in (A). 


Figure 10. Remote sensing signatures during winter and summer. 

Results from co-located active and passive remote sensing instruments (Table 2) viewing 
similar ice and snow conditions (Figure S4). Left panels: measurements during a warming and 
storm event in November 2019; and right panels: during a melting event in September 2020. 
(A) Air temperature and wind speed from the Polarstern weather station and snow surface 
temperature from the IR camera at the remote sensing site (dashed blue line shows time 
periods with potential icing on the lens). (B) Radar backscatter at VV polarization from 
microwave scatterometers L-SCAT at 1.3 GHz and Ku/Ka-radar at 15 and 35 GHz (note the 
different y-scales). (C) Brightness temperature at V polarization from microwave radiometers: 
ELBARA at 1.4 GHz, ARIEL at 1.4 GHz looking at thin ice on a lead, HUTRAD at 7 and 11 GHz, 
Special Sensor Microwave Imager at 19, 37, 89 GHz (not all available data shown). (D) 
Reflected GNSS data, i.e., reflectivity at the remote sensing site (blue) and for sea ice next to 
Polarstern (red). The panel titles give the incidence angles used. Vertical dashed lines mark 
the start of warming and/or storm events. (E) Example photographs of the remote sensing 
site during winter and summer. 


48 


2160 


2165 


2170 


2175 


2180 


2185 


2190 


2195 


2200 


Figure 11. Aerial photo mosaics of the central observatory and its surroundings. 

(A) Thermal infrared image on January 21, 2020, (B) Visual (RGB) image on June 30, 2020. The 
red line shows the perimeter of the central observatory during summer (Leg 4) and its 
approximate position in the winter pack ice during Leg 2 (dashed line), thus allowing the 
alignment of both images. The position of Polarstern at the floe is indicated with the arrow, 
also showing that most parts of the central observatory were at a different location during 
Leg 2 than during Leg 4. The rotation against geographic north resulted mainly from the drift 
of the ice pack past the North Pole, less from rotation. 


Figure 12. Sea ice drift speed as a function of mean wind speed. 

Data were measured onboard Polarstern for one week in January and July 2020, respectively. 
Point clouds show all measurements, large dots and bars give the mean and one standard 
deviation. The free drift estimate (solid line) is based on "Nansen-Ekman ice drift law", which 
assumes that ice drift speed equals to 2% of the wind speed. 


Figure 13. Total deformation derived from two consecutive Sentinel-1 SAR images. 

Data cover 200 x 200 km distance of the MOSAIC central observatory in (A) winter acquired 
on December 30-31, 2019, and (B) summer acquired on June 20-21, 2020. White arrows 
display sea ice motion. 


Figure 14. Sea ice thickness and snow depth distributions. 

Probability density functions (PDF, histogram areas normalized to 1) of (A) total (sea ice plus 
snow) and (B) snow thickness along the transect loops (Figure 2) on January 23 and July 07, 
2020. (C) Melt pond depth on July 07, 2020, as derived from the snow probe modified for melt 
pond measurements. 


Figure 15. Sea ice physical properties from ice cores. 

(A) Stratigraphy of an ice core collected at the first-year ice site on December 02, 2019. Vertical 
(rectangular) and horizontal (circular) thin sections photographed between crossed polarizers. 
(B) Temperature and salinity of first-year sea ice from the ice cores on December 02, 2019 
(solid lines; sample 1_10-FYI, core length 0.71 m), and July 06, 2020, (dashed lines; sample 
4 _46-FYI, core length 1.60 m). Photos of coring sites are shown in Figure S2. The snow/ice 
surface is on top, freeboard at O m. 


Figure 16: Example vertical profiles of the snow pack during winter and summer. 

(A) The snowpack on January 10, 2020 and (B) the sea ice surface with the surface scattering 
layer on July 06, 2020. The depth profiles on the left show density, specific surface area and 
grain size (optically equivalent diameter, OED) analyzed from the micro-CT. In addition, the 3- 
dimensional reconstruction of the microstructure from the same dataset can been seen on 
the right. The y-axis gives the depth with the snow surface at 0 m (top) and the ice surface at 
the bottom of the profile. 
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Overview of the MOSAIC expedition — Snow and Sea Ice 


Table 1. Key dates of the MOSAIC expedition. 


Tables for 


Date Comment Expedition Leg | Observatory ° 
Sep 20, 2019 Departure from Troms@ Start Leg 1 = 

Oct 04, 2019 Start Drift 1 “i Start CO1 

Oct 07, 2019 First buoys deployed in DN1 - Start DN1 
Dec 13, 2019 Resupply and personnel exchange, on site leg1=>Leg2 | - 

Feb 24, 2020 Resupply and personnel exchange, on site Leg2=>Leg3 | - 

May 16, 2020 | End Drift 1 - ; 

Jun 04, 2020 Resupply and personnel exchange, Svalbard Leg3=>Leg4 | - 

Jun 19, 2020 Start Drift 2 - Start CO2 

Jul 31, 2020 End Drift 2 - - 

Aug 12, 2020 Resupply and personnel exchange, Fram Strait | Leg4=>Leg5 | - 

Aug 21, 2020 Start Drift 3 - Start CO3 + DN2 
Sep 20, 2020 End Drift 3 - = 

Oct 12, 2020 Arrival in Bremerhaven End Leg 5 : 


*CO1 to CO3 denote the three different central observatories and DN1 and DN2 denote the two distributed networks. The term “drift” refers to the manned 
drift, when Polarstern was drifting with the corresponding CO. Note that no end dates are given for the COs and DNs, because autonomous stations 
continued reporting beyond the manned drift. 


Table 2. Overview of sea ice remote sensing (RS) instruments operated during MOSAIC. 


Frequency band Name Location Channels Polarization Corresponding satellites/sensors 
Microwave radiometers 
UHF to L-band UWBRAD RS site 0.54, 0.9, 1.38, 1.74 GHz single (right-hand-circular) proposed future mission 
L-band ELBARA RS site 1.4 GHz dual-polarization (H, V) SMOS, SMAP, CIMR 
L-band EMIRAD2 Polarstern 1.4 GHz full polarimetric SMOS, SMAP, CIMR 
L-band ARIEL Transects 1.413 GHz dual-polarization (H, V) SMOS, SMAP, CIMR 
C- to K-band HUTRAD RS site 6.85, 10.65, 18.70 GHz dual-polarization (H, V) AMSR-E/2/3, CIMR 
K- to W-band SSMII RS site 19, 37, 89 GHz dual-polarization (H, V) SSM/I, SSMIS, FY-3, AMSR-E/2/3, CIMR 
Radars 
L-band L-SCAT RS site 1.26 GHz full polarimetric (HH, VV, HV, VH) | ALOS/PALSAR 1, -2, and -3; SAOCOM; NISAR; ROSE-L 
C-band C-SCAT RS site 5.55 GHz full polarimetric (HH, VV, HV, VH) | Sentinel-1, Radarsat-1/2, Radarsat Constellation Mission; ASCAT 
X-band X-SCAT RS site 9.65 GHz full polarimetric (HH, VV, HV, VH) | TerraSAR-X, PAZ, COSMO-Skymed, KOMPSAT-5 
ra iketader———_[ASsteanvamecs | ZJBSHE SOATSHE ag GHIMEET | aronazaon Hv) | SHORE? Seine 3, Aa, ES 17, Eva CHSTAL times 
Others 
; right-handed (RHCP) and left- 
L-band GNSS-R on ice RS site ean eEhaVEeGnz) onde Har ce a ESA FSSCat, UK TDS-1, ESA PRETTY, China FY-3E and Taiwan FS-7R 
(LHCP) 
‘ right-handed (RHCP) and left- 
L-band GNSS-RonPolarstern _| Polarstern So aa Serre Gra) panded ae ee ESA FSSCat, UK TDS-1, ESA PRETTY, China FY-3E and Taiwan FS-7R 
(LHCP) 
Infrared and visual IR/video camera RS site 7.5-14 um, visual Not applicable MODIS, VIIRS, Envisat, Sentinel-3 
Visual to NIR Hyperspectral camera RS site 400-1000 nm in 204 spectral bands Not applicable MODIS, VIIRS, Sentinel-3, Sentinel-2, ICESat-2 


Table 3. Table of key parameters of the MOSAIC snow and sea ice program. 
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o/ 3/2/44); oa;]8;]/sg]0/ 2) S| S|] E 
s/Sla/s/5/z/5/S/S/£/2/ 2/2 
oO/F/2/SiFlose@|asz#/o Sin ie 
Parameter by Category ay Nm) St] OY Mm | Oo] oO) | S| | | Methods or Instruments 
Sea ice and snow mass balance 
Thickness, sea ice - - |X |X ]- - - - - |X |X |X ]- Electromagnetic induction, drilling, 
coring, sonar, thermal properties 
Thickness, snow - |- J- JX |X J- J- J-  |- |X J- |X J]- | Magna probe, laser scanning, drilling, 
coring, sonar, stake readings 
Snow water equivalent - xX |- - - - - - - - - - - | Snow pit 
Bottom topography - - - - - - - - - X |X |- - Electro magnetics, sonar 
Surface roughness - |- J- J- |X Je fe J- fe J-  J- |X J- | Laser scanning 
Freeboard - |- |X JX J- J- f- J-  f- J-  |- |X J- | Coring, laser scanning 
Floe size (distribution) X |- J- J- J- J- J- |X J-  J-  J- |X J- | Cameras, bridge obs. 
Lateral melt rate - |- J- Je Je fe fe fe fe fe f- |X J- | Cameras 
Sea ice and snow physical properties 
Salinity, snow and sea ice - |X |X JX J- J- J- J- J- J-  J-  J-  J- | Snow pit, coring, conductivity sensor, salt 
harp 
Temperature, sea ice - - |X |X j- - - - - - - - - |Temperature probe, thermistor strings 
Density, sea ice - |- |X Je J- fe  f- fe  f- J-  J-  J-  J- | Hydrostatic density kit 
Texture / microstructure, snow and sea ice - - xX |- - - - xX |- - - - - Visual inspection, cameras, Micro-CT 
Porosity, sea ice and ridges - |- |X J-  J- Je f- Je J- |X |X |X J- | Sonar, Electromagnetic induction, coring 
Sea ice microstructure / stratigraphy - X |X |- - - - |X |- - - - - | Coring, Micro-CT 
Snow density, stratigraphy, hardness - |X J- Je Je  J- fe f- fe J-  f-  J-  J- | Snow pit, SMP 
Snow grain size, specific surface area - xX |- - - - - - - - - - - Snow pit, Micro-CT, SMP, NIR 
Temperature, snow - |X J- |X ]- - - - - - - - - |Snow pit, thermistor strings 
Stable oxygen isotope ratio (018), snow and - |X |X J-  J- Je J- J- J-  J-  J-  J-  J- | Snow pit, coring 
sea ice 
Snow surface roughness (small-scale) - |X J- J- Je Je fe fe fe fe Je f-  J- | Snow pit 
Energy budget, optical properties and melt ponds 
Irradiance, spectral and broadband xX |- - - - |X |- - - - |X |X ]- Radiometers (manual, stations, drones, 
helicopter, ROV) 
Albedo, spectral and broadband - - - - - |X |- - - - - |X J- Radiometers (manual, stations, drones, 
helicopter) 
Transmittance, spectral and broadband - - - - - |X |- - - - |X |- - Radiometers (manual, stations, ROV), 
hyperspectral imager 
IOP, spectral, snow and ice and ponds - |- J- J-  J- |X |X J-  J- J-  J-  J-  J- | Optical probe, propelled platform 
Impurities - |X |X ]- - - |X J- - - - - - | Snow pit, coring, water samples 
Depth, melt ponds - |- J- J-  J- J- |X J-  J- |X J- |X ]- | Magna probe water samples, propelled 
platform, cameras 
Geometry and coverage, melt ponds X |- J- J- J- JX J- Je f- J-  J- |X = }- | Cameras 


Temperature, melt ponds 


Water samples 


Salinity, melt ponds 


Water samples 


Bottom imagery, melt ponds 


Cameras 


Surface properties - X |X |Cameras, bridge observations 

Dynamical and mechanical properties 

Location - - |- | Position buoys, GNSS 

Deformation - X |- | Ship Radar, laser scanning, cameras, 
position buoys 

Stress - - |-  |Stress buoys 

Strength - - |- | Uniaxial tests, borehole jack 

Seismics - - - |Seismometers 

Microwave properties 

Emission - - X_ | Radiometers ARIEL, ELBARA, EMIRAD2 
HUTRAD, SSMI, UWBRAD 

Backscatter - - |X |Radars L-, C-, X-SCAT, Ku/Ka-Radar 

Reflectivity - - |X |GNSS-R on ice and Polarstern, Ku/Ka- 
Radar 

Dielectric Permittivity - - |X |Vector Network Analyzer and open-end 
microwave resonator sensor, 
HydraProbe, GNSS-R 

Supplemental observations 

Temperature, ocean X - |-  |Temperature sensors, thermistor strings 

Salinity, ocean - - |- | Conductivity sensor 

Conductive heat fluxes X - |-  |IMB 

Air temperature X - - |Temperature sensors 

Photo documentation - X |X | Cameras 

Infrared images - X |X _ |IR cameras 

General snow and ice conditions - X |- | Bridge observations, cameras 


® Task and method terms are explained in Section 2 with the sub-section given through the number of the task (e.g., “1 General” > “2.1 General snow and 


ice observations”). 


> Entries are sorted by topics according to Section 3 of this manuscript. 
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Figure 1. Drift tracks of the central observatories (CO) of MOSAIC in 2019-2020. 
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Figure 2. Main sites and installations in the central observatory at the end of Leg 3 (CO1). 
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Figure 3. Snow and sea ice observations during the field phase of MOSAIC. 
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Figure 4. Schematics of the snow and sea ice processes studied during MOSAIC. 


Figure 5. The ‘sea ice clock’. 
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Figure 6. Snow redistribution observed via terrestrial laser scanning (TLS). 
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Figure 7. Light transmittance through sea ice measured with the remotely operated vehicle. 
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Figure 8. The Jaridge ridge and adjacent lead in summer 2020 (Leg 4). 
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Figure 9. Surface albedo during summer (Leg 4). 
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Figure 10. Remote sensing signatures during winter and summer. 


10 


Wind speed (m s’) 


A) Jan 21, 2020 B) June 30, 2020 
Winter: Thermal infrared image Summer: True-color image 


GO SEI, Ril SSE 
ea. +; 'y; ie 


ms 


EPSG: 3996 
EPSG: 3996 


M235 | 23602—Ss« 237 238 EE 239 
Temperature (K) 


Figure 11. Aerial photo mosaics of the central observatory and its surroundings. 
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Figure 12. Sea ice drift speed as a function of mean wind speed. 
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Figure 14. Sea ice thickness and snow depth distributions. 
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Figure 15. Sea ice physical properties from ice cores. 
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Figure 16: Exemplary vertical profiles of the sea ice surface and snow during winter and 
summer. 
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S1 Sea ice and snow work and concept 

The realization of the field program by the highly interdisciplinary ICE team required an 
extraordinary amount of coordination and interaction. This effort succeeded and initiated 
manifold new collaborations and scientific exchange. One of the main challenges was to 
ensure consistency in methodology and data quality of the individual observations over the 
year, typically including four to five different principal investigators on board and large teams 
on land, covering a broad range of expertise and specific questions to the same data set. As a 
result, the task structure, as described in Section 2, is not fully consistent as it merges topics, 
instruments and ice types, but was found to be most practical in the organization of the daily 
work in the field. This structure mostly represents how different groups structured their work 
(operational principles) and data sets. 


During the field phase, the ICE team was represented with 12 (Leg 1), 14 (Leg 2), 11 (Leg 3), 
14 (Leg 4), 11 (Leg5) berths on Polarstern and with 9 berths on Akademik Fedorov (Leg 1a). In 
addition, the ICE team was represented by (co-) cruise leaders on Legs 1 and 2 on board of 
Polarstern and Leg 1a on Akademik Fedorov. Overall, 66 different persons participated on both 
ice breakers. Beyond direct participation on board, strong support was given from land before, 
during and after the field experiment. Overall, approx. 150 people were involved in the snow 
and sea ice work and contributed in various ways to the planning, design and successful 
completion of the work program. 


A particular preparation phase was based on the broad experience and expertise across the 
participating researchers. A main process was the agreement on observational protocols for 
all field tasks and methods prior to the expedition. Additional workshops and training 
programs for team members were critical components in preparation for the field experiment. 
Week-long field training courses were held to cross train team members on the full suite of 
snow and sea ice measurement protocols and to perform dedicated instrument tests. These 
courses were held in Hailuoto, Finland (February 22 to March 07, 2019) and Utqiagvik, Alaska 
(April 07 to 13, 2019). More specialized trainings were performed for flight training and system 
testing (particularly for navigation systems) for the unmanned aerial systems near 


Longyearbyen in April 2019. Various cross-calibration initiatives were also realized before and 
after the field phase. 


During the drift, the daily work was organized along weekly plans with designated time slots 
per task. Each team member was assigned specific tasks on each day of the week. The same 
task was mostly carried out by the same people throughout each leg to maintain the highest 
possible consistency. The work on board was supported by the task members at home, 
especially by those who carried out the same measurements on earlier or later legs. Snow and 
sea ice tasks in the field were often supported by volunteers from the other scientific teams 
as well as by the logistics team. Intense observation periods or case studies were incorporated 
into the weekly routine to increase the spatial and/or temporal resolution of existing tasks or 
to capture certain events. A more detailed definition of events during MOSAIC is under 
development by the project coordinators. An example of an intense observation period is also 
discussed in Section 4.2. 


Figure S1 shows how the work from the different tasks was distributed over the respective 
week. In addition to the task work (colored time slots), approximately 1/3 of the time slots 
were used for data work, basic tasks incl. supplemental observations, additional work without 
immediate task relation (FLEX time), and free time. The week of July 06, 2020, included an 
intensive observation period of a 24-hour continuous sampling to capture the diurnal cycle 
under polar day conditions, complementing a similar study during polar night on Leg 1. The 
different schedules for both case studies demonstrate the general realization of the work 
program organized through specific tasks over the entire year. The comparison of these two 
weeks shows how the weekly plans changed over seasons: additional optical measurements 
(‘OPTICS’) were carried out under daylight conditions. The rapid changing surface conditions 
daily snow pit (including surface properties), very frequent ‘TRANSECT’ and additional ‘POND’ 
work in July. The remote sensing work (‘REMOTE’) was allocated more time in January, for 
example when a larger suite of instruments was operated on the ice. Different project and 
process related foci were realized on individual legs like ‘RIDGE’ and ‘POND’ work in summer 
(Leg 4) or ‘DYNAMICS’ work in winter (Leg 2). The concept of ‘FLEX’ time, time that was not 
pre-allocated before the respective leg, turned out to be most beneficial and at the same time 
essential to enable the planned work program. This time allowed reaction on the continuous 
changes and challenges in the field. Also allocating sufficient time for in-field data 
documentation and early processing, as well as dedicated time slots for laboratory work were 
budgeted and needed. Pre-assigned half days off on Sunday mornings, or on other days when 
applicable, paid off given the long field phases and continuous high workloads. 


Week 20-26 January 


Morning 
Monday Tuesday Wednesday Thursday Friday Saturday Sunday 
DATA DATA DATA DATA DATA DATA FREE 
DATA DATA DATA DATA DATA FREE 
DATA DATA DATA DATA DATA FREE 
DATA LAB DATA LAB 


DATA 
HELI HELI 


FREE 
FREE 
FREE 


Afternoon 
Monday Tuesday Wednesday Thursday Friday Saturday Sunday 
DATA DATA FLEX HELI BASIC DATA FLEX 
DATA DATA FLEX DYNAMIC DATA DYNAMIC FLEX 


FLEX DYNAMIC FLEX 
BASIC DYNAMIC FLEX 
DATA FLEX 
DATA FLEX 
DATA FLEX 
DATA FLEX 
DYNAMIC FLEX 
DYNAMIC FLEX 
FLEX 
FLEX 
FLEX 
FLEX 


Week 6-12 July 


Morning 
Monday Tuesday Wednesda Thursday Friday Saturday Sunday 
FLEX FLEX FLEX FLEX FLEX FLEX FREE 
BASIC BASIC BASIC FLEX FLEX FLEX FREE 
BASIC BASIC BASIC FLEX FLEX FLEX FREE 


HELI BASIC FLEX BASIC FREE 
BASIC FLEX BASIC FREE 
LAB BASIC LAB FREE 
HELI BASIC 
DYNAMIC HELI 


DYNAMIC OPTICS 


Afternoon 
Monday Tuesday Wednesday Thursday Friday * Saturday ** Sunday 
FLEX FLEX FLEX FLEX FLEX 
FLEX BASIC FLEX FLEX FLEX 
FLEX DATA FLEX BASIC FLEX 
BASIC HELI BASIC DATA FLEX 


BASIC DATA HELI 
BASIC HELI 

DATA 

LAB 

HELI 


DYNAMIC 
*24-hour sampling period began; **24-hour sampling period ended 


Figure $1: Weekly work plans for snow and sea ice observations. 

Schematics from the weeks starting January 20, 2020, (left) and starting July 6, 2020, (right). 
Each line represents one person. Workdays were split into morning and afternoon blocks. 
Colors are consistent with the sites in Figure 2. Abbreviations refer to the tasks, as described 
in Section 2 of the main manuscript, in addition the following terms are used: ‘FLEX’ for flexible 
tasks, ‘DATA’ for data documentation and processing, ‘BASIC’ for routine work, ‘LAB’ for 


(freezer) laboratory work. 


S2 Methods and field set-up details 


The MOSAiC snow and sea ice program was based on a large number of specialized 
instruments and methods, which are usually referred to in abbreviations or acronyms. In 
addition, the field work concept included terminology (and abbreviations) that are unknown 
to many external readers. All these terms are compiled in Table S1 to ease reading of the 
manuscript. Abbreviations for all remote sensing instruments are given in Table 2 in the main 


text; here only names mentioned in the text are listed. 


Table $1: Instrument names and abbreviations used in the text. 


Short name Full name / description 

ALS Airborne Laser Scanner 

ARIEL Airborne Radiometer in L-band / used on a mobile sled on transects 

co Central Observatories (existence of CO1 to CO3) 

DN Distributed Network (existence of DN1 and DN2) 

EM Electro Magnetic 

EM-Bird Helicopter-towed electro-magnetic sounding instrument 

FY| First Year (sea) Ice 

GEM EM induction sounding instrument 

GNSS (-R) Global Navigation Satellite System (and Reflectrometry) 

HELiPOD Helicopter-towed atmospheric sensor suite 

HUTRAD Helsinki University of Technology Radiometer / Microwave radiometer at 
Remote Sensing site 

IMB (Sea) Ice Mass-balance Buoy 

IR Infrared 

LIDAR Light detection and ranging 

Micro-CT X-ray Micro Computer Tomograph 

MOSAIC Multidisciplinary drifting Observatory for the Study of Arctic Climate 

NIR Near Infrared 

ROV Remotely Operated Vehicle 

RS Remote Sensing 

SAR Synthetic Aperture Radar 

SCAT Scatterometer 

SMP Snow Micro Pen 

SSL Surface Scattering Layer 

SYI Second Year (sea) Ice 

TIR Thermal Infrared 

TLS Terrestrial Laser Scanner 

UHI Underwater Hyperspectral Imager 


Sea ice coring 
Photographs of the conditions at the coring sites from winter and early summer are shown in 
Figure S2. The sea ice data used to create Figure 15, is given in Tables S2 and S3. 


Figure S2: Sea ice coring sites. 

Photographs of the coring site during (A) Leg 3 on March 21, 2020, and (B) Leg 4 on June 22, 
2020. The photo from the spring Leg 3 also shows the shelters that were temporarily set up for 
the coring work. 


Table $2: Sea ice core data: winter 

Salinity (S) and temperature (T) profile for first year (FYI) and second year (SYI) of the 1_10 
coring event on December 02, 2019. The upper (z0) and lower (z1) boundary of the salinity 
section depth, and the depth of temperature measurements (z) are given in m relatively to the 
ice surface. The graph is shown in Figure 16. 


FYI SYI 


0.000 0.050 69) 0.025 -15.3 | 0.000 0.050 0.9] 0.025 -83 
0.050 0.100 5.8) 0.075 -13.0 | 0.050 0.100 O05] 0.125  -7.8 
0.100 0.150 59] 0.175 -12.0 | 0.100 0.150 09] 0.225 -6.5 
0.150 0.200 4.1] 0.275 -8.5 0.150 0.200 1.1] 0325 -6.1 
0.200 0.250 3.5} 0.325 -7.4 0.200 0.250 1.2| 0.375 = -6.2 
0.250 0.310 481] 0.375 -7.1 0.250 0.290 24] 0425  -5.8 
0.310 0.360 481] 0.425 -6.4 0.290 0.350 2.3) 0.545 -5.0 
0.360 0.410 481] 0.475 -5.5 0.350 0.400 3.9) 0.645 = -4.3 
0.410 0455 45] 0.525 -5.3 0.400 0.450 1.7 | 0.745 ~— -3.3 
0.455 0.500 5.5) 0.625 -3.1 0.450 0.500 1.2} 0.845  -1.9 
0.500 0.550 4.6) 0.675 -2.3 0.500 0.550 3.1} 0.870 -1.7 
0.550 0.600 4.7) 0.695 -1.9 0.550 0.600 1.6 


Table $3: Sea ice core data: summer 

Salinity (S) and temperature (T) profile for first year (FYI) and second year (SY!) of the 4_46 
coring event on July 06, 2020. The upper (z0) and lower (z1) boundary of the salinity section 
depth, and the depth of temperature measurements (z) are given in m relatively to the ice 
surface. The graph is shown in Figure 16. 


FYI SY| 


0.910 0.960 3.8 0.890 0.950 0.1 1.625 -0.5 
0.960 1.010 3.7 0.940 1.000 0.2 1.725 -0.9 
1.010 1.060 3.9 0.990 1.050 0.2 1.825 -1.0 
1.060 1.110 3.5 1.040 1.100 0.2 1.925 -1.2 
1.110 1.160 3.5 1.090 1.150 0.2 2.025 -1.4 
1.160 1.210 3.6 1.140 1.210 0.2 2.125 -1.4 
1.210 1.260 3.5 1.190 1.260 0.3 2.210 -1.2 
1.260 1.310 3.1 1.240 1.310 0.2 
1.310 1.360 3.2 1.290 1.360 1.1 
1.360 1.410 2.9 1.340 1.410 1.7 
1.410 1.460 3.2 1.390 1.460 1.0 
1.460 1.510 2.6 1.440 1.510 2.3 
1.510 1.545 1.6 1.490 1.560 2.0 
1.545 1.605 2.2 1.540 1.610 0.2 


Remotely operated vehicle (ROV) 

The ROV has been operated from different sites due to the dynamic icescape: Site ROV 1.0 
was never used for scientific dives, being immediately replaced by ROV 2.0 during Leg 1 (CO1). 
Site ROV 3.0 was used on CO1 during Legs 2 and 3, sites ROV 4.0 and 4.5 were operated on 
CO2 during Leg 4, and site ROV 5.0 was operated on CO3 during Leg 5. Figure S3a shows 
exemplary photographs of the set up (Leg 2 and Leg 5) and maps to illustrate linkages to other 
measurements during spring. 


Figure $3: ROV observations. 

(A) Remotely operated vehicle (ROV) site consisting of the power hub, the surface unit (white 
hut), the tent over the hole (photo March 14, 2020), (B) inside the ROV tent (photo December 
07, 2019), (C) ROV site consisting of the surface unit and the hole without tent (photo August 
25, 2020), (D) dive track (yellow line) and excluded acoustic navigation fixes (red dots) as 
overlay on the airborne laser scanner topography surface map (dive on February 04, 2020), (E) 
schematic overview of the ROV dive range (green circle) during autumn (Leg 5, aerial photo 
from September 06, 2020). Other sites: TLS area in orange ellipse, transects in orange lines, 
snow and surface studies in yellow patches, other installations and sites as small yellow 
squares. 


On-ice remote sensing 

The concept of the on-ice remote sensing measurements was to observe the same snow and 
sea ice surface, or at least the same surface type, with all instruments. In addition, other 
manual measurements were co-located on the same site (Figure 2) and ice conditions. Figure 
S4 shows the arrangement of the individual sensors around the observation site. 


10m 
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WO7-S 


mounted on sleds 


20-50 m 


Figure $4: On-ice remote sensing concept. 

Conceptual layout used for the Remote Sensing Sites on the MOSAIC ice floe. For comparability 
all instruments looked at similar ice and snow. Physical ice and snow properties were sampled 
in the vicinity. Photographs of the Remote Sensing Site are shown in Figure 10. Additional sea 
ice remote sensing observations were performed from Polarstern. 


Sea ice drift forecasts 

A near-real-time drift forecast product for the MOSAIC floe was provided by the Year of Polar 
Prediction (YOPP; Jung et al., 2016) Sea Ice Drift Forecast Experiment (SIDFEx). Several 
operational forecast centers and institutes contributed drift forecasts in near-real-time for 
lead times ranging from days to a year. For MOSAIC a consensus ensemble forecast product, 
based on the different forecast systems, was used. Typically, about five different short-term 
(7-10 days) single-trajectory forecasts, which in particular use recent wind forecasts for 
driving the ice drift, and one seasonal, climatological forecast were merged into a seamless 
ensemble forecast. The consensus forecasts were provided onboard Polarstern through the 
MapViewer system to support decision making, and on land through an online tool 
(https://sidfex.polarprediction.net), for placing orders for satellite imagery. Figure S5 
exemplifies the consensus forecasts product, showing the forecast issued on February 24, 
2020. 


Beyond the drift phase, forecast products from fully coupled models were provided on a daily 
base as well as ensemble predictions of sea ice conditions for the coming months 
(https://nps.edu/web/rasm/predictions). All these sea ice model applications were most 
supportive for the highly complex logistical operations of the supply vessels. Advancing the 
fidelity of different models with the hierarchy will allow for the development of optimized 
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observational networks, as they may be used for sea ice monitoring or advanced field 
programs in the future. 
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Figure S5: The SIDFEx consensus drift forecast for the MOSAIC central observatory (CO1). 

The forecast of the Sea Ice Drift Forecast Experiment —_(SIDFEx, 
https://sidfex.polarprediction.net/) starts on February 24, 2020, at 00:00 UTC, when 
Polarstern was at the given position (grey square). Thin solid lines denote individual (merged) 
forecast ensemble members and the corresponding thick solid line denotes the ensemble mean 
(centroid). Lines are colored by calendar month. Dots denote daily observed positions of 
Polarstern before (black) and after (colored by calendar month) February 24, 2020. 


Technological challenges 
The year-long operations in the central Arctic resulted in specific (technological and 
methodological challenges), in particular with respect to automated systems. 


Challenges were observed in flying the Mavic and Spectra drones close to the North Pole, 
where operators needed to apply the manual mode because the compass reading was not 
correct. However, similar problems were not encountered using the HELiX drone, albeit at 
lower latitudes. Future studies should consider developing and leveraging advanced 
navigation systems such as the D-GPS navigation employed by the DataHawk unmanned aerial 


) 


system deployed to make atmospheric measurements, to avoid the challenges encountered 
by the Mavic and Spectra copters. In addition to challenges posed by the navigation systems, 
other factors, including fast ice drift velocities, fog, and icing conditions resulted in a difficult 
operating environment for drone systems. Despite these challenges, the systems deployed 
combined to provide unique perspectives on broadband and spectral albedo and their 
evolution during the melt and refreezing seasons. 


Many complex instruments, e.g., most of the remote sensing sensors, were designed for 
shorter campaigns. Operating them continuously for a full year led to some instrument 
failures, which only partly could be repaired in the field, and thus led to data gaps or a 
complete stop of measurements for some channels (Figure 3b). This was partly compensated 
by the large suite of complementary measurements, but a larger pool of spare parts and more 
trained personnel could have reduced these downtimes. 
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